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Abstract — We consider the problem of network coding across 
three unicast sessions over a directed acyclic graph. We consider 
a SISO scenario, in the sense that each source and receiver 
is connected to the networls through a single edge. We adapt 
a precoding-based interference alignment technique, originally 
developed for the wireless interference channel, to the network 
setting. We refer to this approach as precoding-based network 
alignment (PBNA). Similarly to the wireless setting, PBNA 
asymptotically achieves half the minimum cut. Different from the 
wireless setting, network topology may introduce dependencies 
between elements of the transfer matrix, which we refer to as 
coupling relations, and can potentially make PBNA infeasible. 
The goal of this paper is to characterize the classes of networks 
for which PBNA is feasible, by identifying and interpreting the 
minimal set of coupling relations introduced by network topology. 
To this end, first, we identify graph-related properties of transfer 
functions, which are essential in identifying the minimal set of 
coupling relations. Then, using two graph-related properties and 
a degree-counting technique, we greatly reduce the set of possible 
coupling relations to just three. Finally, we interpret the three 
coupling relations in terms of network topology and present a 
polynomial-time algorithm to check the feasibility of PBNA. 

Index Terms — network coding, multiple unicasts, interference 
alignment. 



I. Introduction 

UNICAST flows are the dominant form of traffic in most 
wired and wireless networks today. The problem of 
designing and adopting communication protocols, that utilize 
network resources in an efficient way and result in through- 
put close to network capacity, has been an important topic 
of research and speculation. Ever since the development of 
network coding and its success in characterizing the achievable 
throughput for single multicast scenario [1| |2|, there has 
been hope that the framework can be extended to give useful 
answers to network capacity problems (namely inter-session 
network coding), especially the practical setting of networks 
with multiple unicasts. Indeed, there have been limited suc- 
cesses in this domain, such as the derivation of a sufficient 
condition for linear network coding to achieve the maximal 
throughput, in networks with multiple unicast sessions |[3l pi. 
However, scalar or even vector linear network coding pj l6| 
||7] alone has been shown to be inadequate in characterizing 
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Fig. 1. Analogy between a SISO network employing linear network coding 
and a wireless interference channel, each with three unicast sessions (s;, di), 
i = 1,2, 3. Both these systems can be treated as linear transform systems 
and are amenable to interference alignment techniques. 

the limits of inter-session network coding |8 |, which includes 
the network setup with multiple unicasts. Losing the linear 
coding formulation leaves the problem somewhat unstructured, 
and that has stunted the progress in obtaining improved rates 
or even guarantees for a broad class of networks. 

In this paper, we consider the problem of linear network 
coding for networks with multiple unicast sessions, repre- 
sented by a directed acyclic graph (DAG). Similarly to the 
SISO wireless channel, we consider a class of SISO networks, 
where each source and receiver is connected to the network 
through a single edge (see Fig. [T^). The use of linear network 
coding provides a linear transfer function representation for 
a network with respect to its unicast streams [9]. As 
discussed in [3|, the "interference" caused by one unicast 
stream can significantly affect the achievable rates of other 
streams. pT] proceeds to develop a sufficient but highly re- 
strictive condition for "interference-free" transmission to be 
possible in a network with multiple unicast streams. Due to 
the restrictive condition, it is generally very difficult to find 
network codes that result in a interference-free framework. In 
fact, as proved in |10|, the problem of finding linear network 
codes with given rate requirements for multiple unicasts is 
NP-hard. Thus, only sub-optimal and heuristic methods are 
known today, e.g., methods based on Unear optimization | |TT[ 
| [T2| and evolutionary approaches | [T3| . 

As opposed to the interference-free framework |3), we 
present another alternative approach, which allows the exis- 
tence of interferences at each receiver Our idea originates 
from the observation that under the linear network coding 
framework, a SISO network with multiple unicast sessions es- 
sentially mimic the wireless interference channel (Fig.[T]). This 
analogy between networks and wireless interference channel 
enables us to apply the technique of precoding-based interfer- 
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ence alignment, originally introduced by Cadambe and Jafar 
for wireless interference channels, to wired networks. We 
refer to this approach as precoding-based network alignment, 
or PBNA for short. The use of PBNA significantly simplifies 
the code design problem of network coding, as all the encoding 
and decoding operations are predetermined regardless of the 
network topology. Moreover, using this approach, each unicast 
session is capable of achieving a transmission rate close to 
one half of the minimum cut between each sender and its 
corresponding receiver, which is better than what is achieved 
through time-sharing (i.e., the unicast sessions operate one at 
a time in a cyclic fashion). Here, we define the minimum cut 
in the information-theoretic sense | fT5| - this serves as the 
upper bound on the achievable rate (network capacity) for any 
unicast stream. 

As mentioned above, the situation of several unicast sessions 
over a wireline network can be thought of as being equiva- 
lent to supporting the same unicast sessions over a wireless 
interference channel with the same linear transfer functions. 
Network coding across sessions at intermediate nodes emu- 
lates superposition in the wireless channel. Essentially, the 
entire network can be viewed as a channel, albeit a channel 
that is not given by nature, as it is the case in wireless, 
but determined by our routing and coding decisions. This 
has the advantage that it allows us to control the channel. 
However, it also has the disadvantage that it introduces spatial- 
correlation between end-to-end paths that share links among 
themselves. This correlation is not present in wireless channels 
with high probability. This is a very important difference 
between the setup of multiple unicasts vs. traditional wireless 
interference channel; in our problem, there may be depen- 
dencies between elements of the transfer matrix (also called 
transfer functions) introduced by the graph structure, called 
coupling relations, which might make PBNA infeasible in 
some networks [16|. In contrast, the channel gains in wireless 
channel, being inherently independent from each other, make 
precoding-based interference alignment feasible | [T4[ with high 
probability. Another difference is the fact that while channel 
gains in wireless interference channel come from the real or 
complex field, network coding for multiple unicasts has an 
algebraic flavor to it since the network transfer functions and 
packet operations are in some finite field domain. Therefore, 
traditional interference alignment techniques for interference 
channels cannot be directly applied. In this work, we develop 
systematic alignment approaches and algorithms for network 
coding across multiple unicasts, and we characterize their 
feasibility depending on the network structure. Note that the 
interplay of multiple sessions in a general network setup can 
be very complicated with arbitrary correlation between the 
session streams. As such, we focus on a useful special case: 
SISO networks with three unicast sessions; this is the smallest, 
yet non-trivial, instance of the problem and can be used as a 
building block and better understanding for network coding 
across multiple unicasts. 

Main Contributions. The goal of this paper is to fully 
characterize the networks, for which PBNA is feasible, by 
identifying the minimal set of coupling relations that make 
PBNA infeasible, and interpreting these coupling relations in 



terms of network topology. We make the following contribu- 
tions: 

1 ) We show that network topology introduces special prop- 
erties to the transfer function, which we refer to as 
graph-related properties. 

2) Using two graph-related properties and a simple degree- 
counting technique, in conjunction with the results of 
(TT], we identify the minimal set of coupling relations 
existing in graphs, the presence of which makes PBNA 
infeasible. 

3) We interpret the coupling relations in terms of network 
topology. Based on these interpretations, we present a 
polynomial-time algorithm to check the feasibility of 
PBNA. 

The rest of this paper is organized as follows. In Section [ll] 
we review related work. In Section |lllj we show how to apply 
precoding-based interference alignment to the network setting. 



In Section IV we present an overview of our main results 
regarding the feasibility conditions of PBNA. In Section |V] 
we present two graph-related properties of transfer functions, 
which play important roles in identifying the minimal set of 
coupling relations. In Section VI we discuss the feasibility 



conditions of PBNA. In Section VII we provide interpretations 
of the coupling relations and a polynomial-time algorithm to 
check the feasibility of PBNA. Section VIII concludes the 
paper and summarizes future work. 

II. Related Work 

A. Network Coding 

Network coding was first proposed to achieve optimal 
throughput for single multicast scenario [l] |2| |3|, which 
is a special case of intra-session network coding. The rate 
region for this setting can be easily calculated by using linear 
programming formulations [18|. Moreover, the code design 
for this scenario is fairly simple: Either a polynomial-time 
algorithm p9) can be used to achieve the optimal throughput 
in a deterministic manner, or a random network coding scheme 
|20'| can be used to achieve the optimal throughput with high 
probability. 

In contrast, for inter-session network coding, which includes 
the practical case of multiple unicasts, there have been only 
limited progresses. It was observed that there exist networks 
in which network coding significantly outperforms routing 
schemes in terms of transmission rate |4 |. However, there exist 
only approximation methods to characterize the rate region for 
this setting |21|. Moreover, it was proven that finding linear 
network codes for this setting is NP-hard |10|. Therefore, 
only sub-optimal and heuristic methods exist to construct 
linear network code for this setting. For example, Ratnakar 
et al. pT] considered coding pairs of flows using poison- 
antidote butterfly structures and packing a network using these 
butterflies to improve throughput; Draskov et al. |12| further 
proposed a linear programming-based method to find butterfly 
substructures in the network; Ho et al. p2) developed online 
and offline back pressure algorithms for finding approximately 
throughput-optimal network codes within the class of network 
codes restricted to XOR coding between pairs of flows; 
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Effros et al. p3\ described a tiling approach for designing 
network codes for wireless networks with multiple unicast 
sessions on a triangular lattice; Kim et al. |13 | proposed an 
evolutionary approach to construct linear code. Unfortunately, 
most of these approaches don't provide any guarantee in 
terms of performance. Moreover, most of these approaches are 
concerned about finding network codes by jointly considering 
code assignment and network topology at the same time. In 
contrast, our approach is oblivious to network topology in the 
sense that the design of encoding/decoding schemes is isolated 
from network topology, and is predetermined regardless of 
network topology. The isolation of code design from network 
topology greatly simplifies the code design of PBNA. 

B. Interference Alignment 

The original concept of precoding-based interference align- 
ment was first proposed by Cadambe and Jafar |14| to 
achieve the optimal degree of freedom (DoF) for K-user 
wireless interference channel. After that, various approaches 
to interference alignment have been proposed. For example, 
Nazer et al. proposed ergodic interference alignment p4) ; 
Bresler, Parekh and Tse proposed lattice alignment ^251; Jafar 
introduced blind alignment |26| for the scenarios where the 
actual channel coefficient values are entirely unknown to 
the transmitters; Maddah-Ali and Tse proposed retrospective 
interference alignment |27J which exploits only delayed CSIT. 
Interference alignment has been applied to a wide variety 
of scenarios, including K-user wireless interference channel 
p'4|, compound broad channel |28|, cellular networks p9) , 
relay networks [301, and wireless networks supported by a 
wired backbone |31 1. Recently, it was shown that interference 
alignment can be used to achieve exact repair in distributed 
storage systems p2) ||33). 

C. Network Alignment 

The idea of applying interference alignment to network 
coding for multiple unicasts was first introduced by Das et 
al. p4) | [T6( , where the authors proposed an algebraic formu- 
lation of the feasibility conditions for PBNA for three unicast 
sessions. Ramakrishnan et al. observed that the feasibility of 
PBNA depends on network topology, and conjectured that the 
feasibility conditions proposed in |34j can be simplified to just 
two rational functions |35|. Han et al. further confirmed 
this conjecture for the simple case of n = 1. 

The main difference between this paper and prior work is 
that we fully characterize the class of SISO networks for which 
PBNA is feasible, by identifying the minimal set of coupling 
relations, the presence of which invalidates the use of PBNA, 
and interpreting those relations in terms of network topology. 
One should note that there are potentially many multiple 
unicast network setups, where throughputs greater than half the 
mincut per unicast session may be achievable through linear or 
non-linear network coding schemes. Therefore, in general, our 
scheme is not optimal in the sense of achieving the maximal 
throughput. But this "caveat" is not specific to our approach 
alone. 



Our work is closely related to identifying algebraic proper- 
ties of transfer functions which are closely related to network 
topology. Some of recent work also focused on this topic: 
Ebrahimi and Fragouli p6) found that the structure of a 
network polynomial, which is the product of the determinants 
of all transfer matrices, can be described in terms of certain 
subgraph structures; Weifei et al. p7) proposed the Edge- 
Reduction Lemma which makes connections between cut sets 
and the row and column spans of the transfer matrices. 

III. Applying Precoding-Based Interference 
Alignment to Networks 

The basic idea of network alignment is that under linear 
network coding, the network behaves like a linear system, 
similar to the wireless interference channel.' |34 | This analogy 
enables us to borrow some techniques, such as interference 
alignment f\A\, which were originally developed for the 
wireless interference channel, to the network setting. In this 
section, we first illustrate this analogy by revisiting some basic 
concepts of linear network coding. Next, we show how to 
apply precoding-based interference alignment in the context 
of networks by revisiting the aUgnment scheme proposed by 
Cadambe and Jafar |14|. Then, we formulate the feasibility 
conditions of PBNA, in terms of transfer functions; this forms 
the basis for our further discussion. Finally, we discuss some 
important questions regarding the feasibility conditions of 
PBNA. 

A. Treating Networks as Wireless Interference Channel 

In order to provide a concise and clear description of our 
ideas, we consider the simplest non-trivial SISO network 
scenario of three unicast sessions, where each sender and 
receiver is connected to the network via an edge. We assume 
the network is represented by a delay-free directed acyclic 
graph G — ( V, £'), where V is the set of nodes and E the set of 
edges. Without loss of generality, each edge has unit capacity, 
i.e., can transmit one symbol of finite field Fq in a unit time, 
and represents an error-free channel. For the ith unicast session 
a — 1,2,3), let Si and denote its sender and receiver 
respectively. For the ith unicast session, let Ci and denote 
the edges which link the Si and di to the network respectively. 
Gi and Ti are referred to as the sender edge and receiver edge 
for the ith unicast session. For e E E, let head{e) and tail{e) 
denote the head and tail of e respectively. 

Under linear network coding |2| |3|, the symbol transmitted 
along edge e, denoted by Y^, is a linear combination of 
incoming symbols at tail{e), 

Ulhead{e')=ta^Ke) ^e'eYe' Otherwise. 

where Xi represents the source symbol transmitted at (Ji\ Xe'e 
is a variable, which takes values from and represents the 
coding coefficient used to combine the incoming symbol along 
e' into the symbol along e. Due to the linear operations at each 
edge, the output at is a mixture of source symbols 

= TOh(x)Xi + ?7l2i(x)X2 + ?7l3i(x)X3 (1) 
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Fig. 2. Applying precoding-based interference alignment to the network 
setting: At eacli sender edge Ci (i = 1,2,3), tlie input vector Xi is first 
encoded into 271 + 1 symbols thi'ough the precoding matrix V^; then the 
encoded symbols are transmitted through the network in 2ra + 1 time slots; 
at each receiver edge t;, the undesired symbols are aligned into a single 
linear space, which is hnearly indepdent from the linear space spanned by the 
desired signals, such that the receiver can decode all the desired symbols. 



where the hnear coefficient mji{x.) (j = 1,2,3) is called the 
transfer function from aj to t^. According to |[3], m,ji{x.) is a 
multivariate polynomial in terms of all the coding variables 
X = {xee' '■ e, e' S E',head{e) = tail{e')), with each 
monomial corresponding to a path from aj to r^: 



(2) 



where Vji denotes the set of paths from aj to Ti, and tp{x.) 
is the product of all the coding variables along path P and 
represents a distinct monomial in ■mji{x.). 

Note that Eq. ([T]) is analogous to that of a wireless interfer- 
ence channeQ which is shown below: 



= HuWi + H2^W2 + H3,W3 + N, 



1,2,3 (3) 



where Wj, Hji, Ui, and Ni (j = 1,2,3) are all complex 
numbers, representing the transmitted signal at sender j, the 
channel gain from sender j to receiver i, the received signal at 
receiver j, and the noise term respectively. Apparently, Xj's 
(j 7^ *) pl^y the roles of interfering signals, and transfer 
functions play the roles of channel gains. Hence, we can treat 
the network as a mimic "wireless interference channel," and 
utilize interference alignment in the network. Yet, despite this 
formal similarity, there are two major differences between 
the two settings. First, in networks, the transmitted signals 
belong to a finite field F^, whereas in wireless channel, the 
transmitted signals usually belong to the field of real numbers 
M or complex numbers C, both of which contain an infinite 
number of elements. Moreover, in networks, the coefficients 
of received signals are introduced by the linear operations in 
the middle of the network, and are polynomials in terms of the 
coding variables in the network, whereas in wireless channel, 
the channel gains are introduced by nature, and are inherently 
structureless. Indeed, as we will see later, these two differences 
greatly affect the feasibility of PBNA. 

B. Precoding-Based Network Alignment 

Now we explain how to apply precoding-based interference 
alignment to the above network setting by revisiting the align- 

'The wireless interference channel that we consider here has only one sub- 
channel. 



ment scheme proposed by Cadambe and Jafar p4) . We assume 
all TOij(x)'s — 1,2,3) are non-zeros. The case where 
some TOij(x)'s ii ^ j) are zeros can be dealt with similarly, 
and is deferred to Section VI We use Fig. |2] to illustrate the 
basic idea of PBNA. As shown in this figure, at each sender 
edge CTi {i = 1,2,3), we first group the source symbols to 



■ X. 



ki(n)\T 



where 



be transmitted into a vector = {Xl , ■ 
n is a positive integer, and ki{n) equals n + 1 for i = 1, 
and n for i — 2,3. We then use precoding matrix to 
encode into 2n + 1 symbols, which are then transmitted 
through the network in 2?! + 1 time slots. It is easy to see 
that Vi is a (2n + 1) x (?i + 1) mati'ix, and V2, V3 are both 
{2n + 1) X n matrices. As a result, at each receiver edge Ti, 
the output vector, denoted by Z^, is a combination of input 
vectors: 

Z, = MnViXi + M2,V2X2 + M3,;V3X3 

where Mj^ (j = 1,2, 3) is a diagonal matrix, with the 
{k, k) element being the transfer function mji{x'^), where x*^ 
represents the vector of the coding variables for the kth time 
slot. In general, V; may still contain variables. Let ^ denote 
the vector consisting of all the variables in x^, • • • , x^"+^ and 
Vi, V2, V3. We require the following conditions are satisfied 
under some assignment of values to ^ |14J: 

£/i : span(M2iV2) = span(]V[3iV3) 
J2/2 : span(M32V3) C span(Mi2Vi) 
£/3 : span(M23V2) C span(]V[i3Vi) 
^1 : rank(MiiVi M21V2) = 2n + 1 
^2 : rank(Mi2Vi M22V2) 2n + 1 
^3 : rank(Mi3Vi M33V3) = 2n + 1 

In the above conditions, for a matrix E, we use span(E) 
to denote the linear space spanned by the column vectors 
contained in E. .2/1,^2/2,^3 are the alignment conditions and 
SSi, are the rank conditions. Basically, each alignment 
condition guarantees that the undesired symbols or interfer- 
ences at each receiver are mapped into a single linear space, 
such that the dimension of received symbols or the number 
of unknowns is decreased. Moreover, each rank condition 
guarantees that the linear space spanned by the interferences is 
linearly independent from that spanned by the desired symbols, 
and thus each receiver can decode the desired symbols from 
the received symbols. In summary, when all these conditions 
are satisfied, the three unicast sessions can achieve a rate tuple 



^ A 



n+l 



1^^), which approaches ( 



1 1 1 N 



as 



■^n ~ ^ 2n+l ' 2n+l ' 2n+l ^' V2'2'2' 

n -+ 00. In this case, we say that R* is feasible through 
PBNA. 

Note that the alignment conditions and rank conditions are 
all defined on the field of rational functions Fq(^). This poses 
no problem for the alignment conditions, because whenever 
they are satisfied on Fg(^), they are automatically satisfied 
under any assignment of values to ^. However, this is not 
always true for the rank conditions, unless the size of F^ 
is sufficiently large. In contrast, in wireless channel, the 
alignment and rank conditions are all defined on or 
C(^'), where ^' consists of the channel gains for the 2ri + 1 
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time slots. Due to the infiniteness of M and C, as long as 
the rank conditions are satisfied on or C(^'), they are 

satisfied almost surely under a random assignment of values to 



We further define the following functions: 



^' in the context of wireless channel 1 14 1. This indicates that 



the size of Fg does affect the feasibility of PBNA. To avoid this 
issue, we assume that the size of ¥q is sufficiently large. Let 
(f>ii£,) £ denote the determinant of the matrix involved 

in Due to the Shwartz-Zippel Lemma p8) , if the rank 
conditions are satisfied on Fg(^), i.e., 0(^) = 4>i{(,)(f>2{£,)4'3{0 
is a non-zero polynomial, we can find an assignment of values 
to with high probability such that 0(^) is non-zero under this 
assignment. 

C. Algebraic Formulation of the Feasibility Conditions 

As shown previously, the feasibility conditions of PBNA, 
namely the alignment conditions and the rank conditions, are 
all expressed in matrix forms, which are difficult to analyze 
and reveal little about their connection with the network 
topology. For these reasons, before proceeding, we reformulate 
the feasibility conditions in terms of transfer functions. To this 
end, we first reformulate the alignment conditions as follows: 

ia// : M21V2 = M31V3A 
: M32V3 = M12V1B 
s^l : M23V2 = M13V1C 

where A is an 71 x n invertible matrix, and B, C are both 
(n + 1) X 71 matrices with rank n. A direct consequence of ^2/3 
and is that the precoding matrices are not independent from 
each other: Both V2 and V3 are determined by Vi through 
the following equations: 



V2 =Mi3M23lViC V3 = M12M32IV1B 



(4) 



Substituting the above equations into the three alignment 
conditions can be further condensed into a single equation: 

TViC = ViBA (5) 

where T = Mi3M2iM32Mj"2^M^3^M3j^^. Eq. ^ suggests 
that alignment conditions introduce self-constraint on Vi. 
Thus, in general, we cannot choose Vi freely. Indeed, Eq. (jSj) 
is also the major restriction on Vi. Finally, using Eq. Q and 
Eq. (|5]l, the rank conditions are transformed into the following 
equivalent equations: 

^[ : rank(Vi PiViC) =2n + l 
SS'^ : rank(Vi P2V1C) = 2n + 1 
^'3 : rank(Vi P3V1CA-I) = 2n + I 

where Pi = MigMziMlfi^M^g^, P2 = Mi3M22Mj"2^M^3\ 
and P3 = M2iM33M^3^M3i^ Thus, we have simplified 
the feasibility conditions into four equations, i.e., Eq. (jSj) 
and ,^i,^2j=^3- Note that all these equations are formulated 
in terms of Vi and Pj. We summarize this result into the 
following lemma: 

Lemma III.l. Assume all mij{x.)'s = l,2,3j are non- 
zeros. R* is feasible through PBNA if and only if 1) Eq. ([5]) 
is satisfied, and 2) i^'^T^^T^'z '^^^ satisfied. 



Pi(x) = 
P3(x) = 



mi3(x)m2i(x) 

'Till (x) 77123 (x) 
'•»2l(x)77l33(x) 
77723 (x)77l3l(x) 



P2(x) 

?7(x) = 



77Zi3(x)7W22(x) 

TOl2(x)7n23(x) 
77713 (x)77l2l(x)77732(x) 
77712 (x)7n23(x)?7l3i(x) 



(6) 



Clearly, Pi(x) and 77(x) form the elements along the diagonals 
of Vi and T respectively. 

Next, we reformulate the feasibility conditions in terms of 
Pi(x) and 7;(x). To this end, we need to know the internal 
structure of Vi. We distinguish the following two cases: 

Case I: r]{x) is non-constant, and thus T is not an identity 
matrix. For this case, Eq. Q becomes non-trivial, and we 
cannot choose Vi freely. We use the following precoding 
matrices, which were first proposed by Cadambe and Jafar 
03: 

(w Tw • • • T"w) (7) 
= Mi3M23^(w Tw • • • T"~iw) (8) 
V3 = Mi2M32^(Tw • • • T"w) (9) 

where w is a column vector of 2?7 + 1 ones. It is straightfor- 
ward to verify that the above precoding matrices satisfy the 
alignment conditions. Note that the above precoding matrices 
correspond to the configuration where A = I„, C consists of 
the left 77 columns of I„+i, and B the right n columns of 

In order to reformulate the rank conditions, we consider the 
following matrix. 



H 



//i(y^) /2(y^) ••• friy')\ 
/i(y') /2(y') ••• fAy^) 



\fi{yn /2(y0 



friynJ 



where fi{y) (i — 1, 2, • • • , r) is a rational function in 
¥q{y), and the jth row of H is simply a repetition of the 
vector (/i(y); • • • , /r (y)), with y being replaced by y-'. Due 
to the particular structure of H, the problem of checking 
whether H is full rank can be simplified to checking whether 
/i(y)i ■ ■ ' 7 fr{y) are linearly independent, as stated in the fol- 
lowing lemma. Here, /i(y), • • • , /r(y) are said to be linearly 
independent, if for any scalars ai, • • • ,arG Vg, which are not 
all zeros, ai/i(y) H h arfr{y) 7^ 0. 

Lemma III.2. dct(H) ^ Q if and only if fi{y), - • ■ , fr{y) 
are linearly independent. 

Proof: See Theorem 1 of p7) . ■ 
An important observation is that using the precoding ma- 
trices defined in Eq. (|7|-(j9]l, all of the matrices involved in 
^'i, ^3 have the same form as H. Specifically, each row 
of the matrix in 



(1 77(X) 



is of the form: 
•• 7;"(x)p,(x) ••• p,(x)77"-i(x)) (10) 



Hence, using Lemma III.2[ we can quickly derive: 



Lemma III.3. Assume all 'mij{x.)'s (i,j — l,2,3j are non- 
zeros, and 77 (x) is non-constant. R* is feasible through PBNA 
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if for each i — 1,2, sj^ 

P^(x) iSn^ ■ /W'^W e ¥,[z]J{z)g{z) + 0, 



icd(/(z),.g(z)) = < n, 



dg <n 



(11) 



Proof. If Eq. ( 1 1 1 is satisfied, the rational functions in Eq. 
( [Tol l ^re linearly independent. Therefore, due to Lemma |III.2| 
condition is satisfied. Meanwhile, the precoding matrix 
Vi defined in Eq. (j?]) satisfies Eq. (|5]l. Hence, R* is feasible 

■ 

/(')(x)) 



through PBNA by Lemma III.l 



/('?(x)) 



S Sn represents 
the violation of 



Note that each rational function 

a constraint on p,{x), i.e., pi{x) ^ ^(^^, 
which invalidates the use of the PBNA through the precoding 
matrices defined in Eq. (|7]l-(|9]). Also note that Eq. ( [TT] i only 
guarantees that PBNA is feasible for a fixed value of n, i.e., 
each unicast session only achieves a transmission rate close to 
one half of its minimum cut. In order for each unicast session 
to asymptotically achieve a transmission rate of one half of its 
minimum cut, we simply combine the conditions of Lemma 



III. 3 for all possible values of n, and get the following result: 



Theorem III.l. Assume all mij{-x)'s = l,2,3j are non- 
zeros, and ri{x.) is non-constant. The three unicast sessions can 
asymptotically achieve the rate tuple ^, ^) through PBNA 
if for each i — 1,2,3, 



p,(x) i S' 



[ /(^(x)) 



/(z),.g(z)eF,[z],/(z).g(z)^0, 



icd(/(z),.g(z)) = l 



Proof: Clearly, if Eq. (12i is satisfied 



(12) 

R* is feasible 

through PBNA for all possible values of n. Thus, each unicast 
session can asymptotically achieve one half rate as n — > oo. 

■ 

Case II: rj{x.) is constant, and thus T is an identity matrix. 
For this case, Eq. ^\ becomes trivial. In fact, we set BA = 
C, and hence Eq. d5]l can be satisfied by any arbitrary Vi. 
Specifically, we use the following precoding matrices: 



Vl — (%)(2n+l)x(n+l) 
V2 - Mi3M23^ViI' 
V3 = M12M32IV1I' 



(13) 

(14) 

(15) 



where 9ij is a variable, taking values from Fg, and I' consists 
of the left n columns of I„+i. Note that the above precoding 
matrices correspond to the configuration where A = I„, and 
B = C = I'. Clearly, the above precoding matrices satisfy 
the alignment conditions. Meanwhile, using these precoding 
matrices, each row of the matrix in 3§l is of the following 
form: 



'1 C2 



Pi(x)6'i • • • _Pi(x)0„) 



(16) 



Hence, using Lemma III. 2 we can quickly derive 



-Notation: For two polynomials f{x) and g{x), let gcd(f{x),g(x)) denote 
their greatest common divisor, and df the degree of f{x). 




(a) pi(x) : 
r?(x) = 1) 



P2(x) = P3(x) 



(b) pi(x) 



l + r,(x) 



Fig. 3. Examples of realizable coupling relations: The left network realizes 
the c oupling relations Pi(x) = »7(x) = 1 such that the conditions of Theorem 



III.2 



are violated; in the right network, »y(x) 7^ 1, but pi(x) 



Which violates the conditions of Theorem IIII. 1 



l + r;{x)' 



Lemma III.4. Assume all mij{x)'s (i,j ~ l,2,3j are non- 
zeros, and rj{-x) is constant. R* is feasible through PBNA if 
for each i ~ 1,2, 3, Pi(x) is not constant. 



Proof Apparently, the functions in Eq. (I61 are linearly 
independent, and therefore is satisfied due to Lemma 
III. 2 Moreover, Eq. (|5]l is also satisfied. Thus, R* is feasible 



through PBNA by Lemma III.l 



Similarly, the above theorem can be easily extended to the 
asymptotic case: 

Theorem III.2. Assume all jnij(x)'s — 1,2, 3j are non- 
zeros, and ?7(x) is constant. The three unicast sessions can 
asymptotically achieve the rate tuple (i, ^, ^) through PBNA 
if for each i = 1, 2,3, Pi(x) is not constant. 



Proof: By Theorem III.2 if Pi(x) is non-constant, R* is 



feasible through PBNA for any value of n. Thus, each unicast 
session can asymptotically achieve one half rate as n — > 00. 



D. Further Discussion 

In the previous section, we reformulated the feasibility 
conditions in terms of the transfer functions Pi(x) and 7?(x). 
The critical question is: what is the connection between the 
reformulated feasibility conditions and network topology? 

We start by illustrating that through examples of networks 
whose structure violates the feasibility conditions. Let's first 
consider the network shown in Fig. [3a] Due to the bottleneck e, 
it can be easily verified that pi (x) =p2(x) =p3(x) = 77(x) = 
1, and thus the conditions of Theorem |III. 2 [ are violated. More- 
over, consider the network shown in Fig. 3b It can be easily 



verified that for this network, 77(x) 7^ 1, and pi(x) = 
Thus the conditions of Theorem |III.1| are violated. Moreover, 
by exchanging cti o CT2 and ti o T2, we obtain another 
example, where P2(x) = 1 -\-r]{yL), and thus the conditions of 
Theorem [III. 1 1 are again violated. While the key feature of the 
first example can be easily identified, it is not obvious what 
are the defining features of the second example. Nevertheless, 
both examples demonstrate an important difference between 
networks and wireless interference channel: In networks, due 
to the internal structure of transfer functions, network topology 
might introduce dependence between different transfer func- 



tions, e.g., Pi(x) = 1 or Pi(x) = 3^^^^^; in contrast, in 
wireless channel, channel gains are inherently structureless, 
and thus can change independently . 
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Note that the above dependence relations between transfer 
functions can be generalized to the following form: 



(17) 



where f{zi, Z2, - ■ ■ , Zk) is a polynomial in Fg[zi, • • • , Zk]- We 
call s uch relation a coupling relation. As shown in Theorem 



III.l 



G S' represents a coupling 
Given a coupling relation, if there 



each rational function 

relation pifx) = ^^7^4^. 
are networks for which it holds, we say that it is realizable. 

The existence of coupling relations greatly complicates the 
feasibility problem of PBNA. As shown previously, most of 
the coupling relations, such as pi(x) = 1 andpi(x) — j^if^—^, 
are harmful to PBNA, because their presence violates the 
feasibility conditions. The only exception is 77(x) = 1, which 
does help simplify the construction of precoding matrices, 
and thus is beneficial to PBNA. Indeed, as we will see 



in Section IV this coupling relation allows interferences to 



be perfectly aligned at each receiver, and thus each unicast 
session can achieve one half rate in exactly two time slots. 



Unfortunately, as we will see in Section VII this coupling 
relation requires that the network possesses peculiar structure, 
and is generally unavailable in most networks. For this reason. 



we will mainly focus on Theorem III.l which is applicable 
for most networks. This is in stark contrast with wireless chan- 
nel, where precoding-based interference alignment is feasible 
almost surely due to the structureless nature of channel gains 

HE)- 

One interesting observation is that not all coupling relations 
are realizable. For example, consider the coupling relation 
Pi(x) = ?y^(x), where both pi(x) and r/(x) are non-constants. 
Let pi(x) = ^y^, ''7(x) = denote the unique forms^of 
Pi(x) and 77(x) respectively. Consider a coding variable Xee' 
that appears in both and Because the maximum 

degree of each coding variable in a transfer function is at most 
one, according to Eq. (j6|, the maximum of the degrees of Xee' 
in u(x) and t;(x) is at most two. However, it can be easily seen 
that the maximum of the degrees of x^e' in s"^ (x) and t'^ (x) is 
at least three. Therefore, it is impossible that pi(x) — ifix). 
This example suggests that there exists significant redundancy 



44^ e S' are 



in the feasibility conditions of Theorem III. 1 More formally, 
it raises the following important question: 
Ql: Which coupling relations Pi(x) = 
realizable? 

The answer to this question allows us to reduce the set S' 
defined in Theorem |III. 1 1 to its minimal size. For i = 1,2,3, 
we define the following set, which represents the minimal set 
of coupling relations we need to consider: 



is realizable > (18) 



Then the next important question is: 

Q2: Given Pi(x) = G S^, what are the defining 

features of the networks for which this coupling relation holds? 

As we will see in the rest of this paper, the answers to Ql 
and Q2 both lie in a deeper understanding of the properties of 

'For a non-zero rational function h{y) £ ^qiy), its unique form is defined 
as h{y) = where f{y),g{y) e F,[y] and gcd(/(y), g(y)) = 1. 



transfer functions. Intuitively, because each transfer function 
is defined on a graph, it usually possesses special properties, 
called graph-related properties, which are closely related to 
the graph structure. The graph-related properties not only 
allow us to reduce S' to the minimal set S'^, but also enable us 
to identify the defining features of the networks which realize 
the coupling relations represented by 

In the above discussion, we have overlooked a subtle 
issue. In the derivation of Theorem |III.1[ we only consider 
the precoding matrices defined in Eq. (|7]l-(j9|l. However, the 
choices of precoding matrices are not limited to these matrices. 



In fact, as we will see in Section VI given different A, B, and 
C, we can derive different precoding matrix Vi such that Eq. 
(|5]l is satisfied. This raises the following interesting question: 
Q3: Assume some coupling relation Pi(x) = G SI 

is present in the network. Is it still possible to utilize PBNA 
via other precoding matrices instead of those defined in Eq. 

0W 



As we will see in Section VI the answer to this question 



is negative. The basic idea is that each precoding matrix Vi 
that satisfies Eq. dS]) is closely related to the precoding matrix 
defined in Eq. (|7]i through a transform equation. Using this 
transform equation, we can prove that if the precoding matrices 
cannot be used due to the presence of a coupling relation, then 
any precoding matrices cannot be used. 

IV. Overview of Main Results 
In this section, we state our main results regarding the 



feasibility of PBNA. Proofs are deferred to Section VI and 
Appendix|B] Since the construction of Vi depends on whether 
?7(x) is constant, we distinguish two cases. 

A. ?7(x) Is Not Constant 

Theorem IV.l (The Main Theorem). Assume all mij{x)'s 
(i,j = 1,2,3) are non-zeros, and ?7(x) is not constant. The 
three unicast sessions can asymptotically achieve the rate tuple 
|, through PBNA if and only if the following conditions 
are satisfied: 



Pi(x)^5;= l,r,(x), 



?7(x) 



1 + 77(x) 
P2(x)^5^ = {l,7?(x),l + r;(x)} 
P3(x)^5^ = {l,7?(x),l + r;(x)} 



(19) 

(20) 
(21) 



Note that in the Main Theorem, we reduce the feasibility 
conditions of Theorem |III. 1 1 to its minimal size, such that each 
5j' (i = 1, 2, 3) represents the minimal set of coupling relations 
that are realizable. Moreover, as we will see in Section |VII| 
each of these coupling relations has a unique interpretation 
in terms of the network topology. The interpretations of these 
coupling relations further provide a polynomial-time algorithm 
to verify the feasibility of PBNA. 

Remark: The conditions of the Main Theorem can be 
understood from the perspective of interference channel. As 
shown previously, under linear network coding, the network 
behaves as a 3-user wireless interference channel, where the 
channel coefficients mij{x) are all non-zeros. Let H denote 



8 



the matrix with the (i, j)-element being m.ij{x.). Suppose that 
there exists a (fcJ)-Minor of H (fc 7^ which equals zero, 
i.e., A/j./(H) = 0. For such a channel, it is known that the 
sum-rate achieved by the three unicast sessions cannot be 
more than 1 in the information theoretical sense (see Lemma 
1 of |39J), i.e., no precoding-based scheme, linear or non- 
linea:]^ can achieve more than 1/3 per user Therefore, given 
that channel coefficients are always non-zero, the condition 
that Mki(H.) 7^ for k I is information theoretically 
necessary for rate 1/2 feasibiUty. Note that the conditions 
pi(x) 7^ l,pi(x) 7^ 77(x),pi(x) ^ can be respectively 

rewritten as follows: 



mii(x) 7^ 
mii(x) 7^ 
mii(x) 7^ 



mi3(x)m2i(x) 

"123 (x) 
mi2(x)m3i(x) 

m32(x) 

mi3(x)rn2i(x) 



mi2(x)m3i(x) 
"132 (x) 



"123 (x) 

The other conditions can be reformulated by permuting the in- 
dices. Apparently, the first two are equivalent to M32(H) 7^ 
and il/23(H) 7^ respectively. Hence, they are simply the 
information theoretic necessary conditions, so they must hold 
for any achievable precoding-based scheme. The third condi- 
tion involves a sum of the first two, which is not necessary 
information theoretically, but is needed for the precoding- 
based scheme we utilize. In fact, consider the network shown 



in Fig. 3b where the third condition is violated. It is fairly easy 



to find a linear network code such that 1/2 rate is achieved. 

B. ry(x) Is Constant 

In this case, we can choose Vi freely by setting BA = C, 
and thus the feasibility conditions of PBNA are significantly 
simplified. Moreover, each unicast session can achieve one 
half rate in exactly two time slots, as stated in the following 
theorem: 

Theorem IV.2. Assume all mij(jK.)'s (i,j = l,2,3j are non- 
zeros, and ?7(x) is constant. The three unicast sessions can 
achieve the rate tuple ^, |) in exactly two time slots 
through PBNA if and only if for each i — 1, 2, 3, Pi(x) 7^ 1. 

V. Graph-Related Properties 

Since the transfer functions are defined on graphs, they 
exhibit special graph-related properties introduced by the 
graph structure. In this section, we present two such properties, 
namely Linearization Property and Square-Term Property, 
which play central roles in the proof of the Main Theorem. We 
defer the proofs of these two properties to Appendix |A] The 
other graph-related properties are related to the interpretation 
of coupling relations defined in Eq. ([T9]l-(21 1 and are deferred 
to Appendix |A[ 

"^In the precoding-based scheme mentioned here, the encoding/decoding 
process at each sender/receiver might be Unear or non-hnear, but the opera- 
tions at each internal node are linear. 



In the following discussion, we consider the general form 
of the transfer function Pi{x), as defined below 

mQb(x)TOpq(x) 



(22) 



ma5(x)mpb(x) 

where a,b,p,q = 1, 2, 3 and a ^ p,b ^ q. Moreover, by the 
definition of transfer function, the numerator and denominator 
of h{x) can be expanded respectively as follows: 

m,,(x)m,,(x) = E(p,.p,)ep„,,p^, ip.(x)ip.(x) 



iP3(x)ip,(x) 



Hence, each path pair in Tab x Vpq contributes a term in 
ma6(x)mpg(x), and each path pair in Vaq x Vpb contributes 
a term in maq{:x-)mpb{x.) . 

A. Linearization Property 

The first graph-related property, namely Linearization Prop- 
erty, is stated in the following lemma. According to this 
property, if Pi{x) 7^ 1, it can be transformed into its simplest 
non-trivial form, i.e., a linear function or the inverse of a linear 
function, through a partial assignment of values to x. 

Lemma V.l (Linearization Property). Assume ft,(x) is not 
constant. Let h{x) — such that gcd(M(x), w(x)) = 1. 

Then, we can assign values to x other than a variable 
Xee' such that ■u(x) and vipc) are transformed into either 

u(Xee') = CiXee'+Co, v{Xee') = Ci Or u{Xee') = C2,v{Xee') = 

ciXee' +Co, where cq, Ci, C2 are constants in ¥q, and C\Ci 7^ 0. 

The key to the above lemma is to find a subgraph H and 
consider /i(x) restricted to H, i.e., /i(xh) = "'^'>(^h);"p.(^^) 

^ ' y y \ n/ m„„ (xh jnipb (xff ) ' 

where consists of the coding variables in H. In fact, due to 
graph structure, we can always find H such that some variable 
Xee' appears exclusively in the numerator or the denominator 
of h{xH)- Thus, by assigning values to xh other than Xee', 
we can transform h{xH) into a linear function or the inverse 
of a linear function in terms of Xee'- Since h{xff) can be 
acquired through a partial assignment to x, this transformation 
also holds for the extended graph G' . 



B. Square-Term Property 

The second property, namely Square-Term Property, is pre- 
sented in the following lemma. According to this property, the 
coefficient of x^^, in the numerator of /i(x) equals its counter- 
part in the denominator of /i(x). Thus, if x^g, appears in the 
numerator of h{x} under some assignment to x, it must also 
appear in the denominator of h{x), and vice versa. 

Lemma V.2 (Square-Term Property). Given a coding vari- 
able Xee', let /i(x) and /2(x) be the coefficients of x\^i 
in mah[y^'rfipq(x) and maq{'X-)mpb{x) respectively. Then 

/l(x)=/2(x). 

The key idea of the above lemma is that due to graph 
structure, each path pair which contributes an x^^, term in 
the numerator of h{x) corresponds to another path pair which 
contributes an equivalent x^^, term in the denominator of h{x). 
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This correspondence relation automatically yields a one-to- 
one mapping from the x^^, terms in the numerator of h{x) to 
those in the denominator of h{x.). Thus, the summation of the 
Xgg, terms in the numerator of /i(x) equals the summation 
of the Xgg, terms in the denominator of /i(x), and hence 

/l(x)=/2(x). 

VI. Feasibility Conditions of PBNA 
In this section, we explain the main idea behind the proof of 



our main results presented in IV Consistently with Section IV 



we distinguish two cases based on whether 77 (x) is constant. 

A. 77 (x) Is Not Constant 

In this subsection, we first present a simple method to 
quickly identify a class of networks, for which PBNA is 
feasible. Then, we explain how to reduce the set S' defined 
to the minimal set S' 



in Theorem 



III.l 



defined in the Main 
Theorem. This answers the first question we raise in Section 
III-D and automatically proves the "if" part of the Main 



Theorem. Next, we provide an answer to the third question 



we raise in Section III-D which actually proves the "only if" 
part of the Main Theorem. 

As shown 



1) A Simple Method Based on Theorem III.l 



in Theorem III.l the set S' contains an infinite number of 



rational functions, and thus it is impossible to check the fea- 



sibility conditions of Theorem III. 1 in practice. Interestingly, 
the theorem directly yields a simple method to quickly identify 
a class of networks for which PBNA is feasible. The major 
idea of the method is to exploit the asymmetry between pi (x) 
and ri{x.) in terms of effective variables. Here, given a rational 
function /(y), we define a variable as an effective variable of 
/(y) if it appears in the unique form of /(y). Let V(/(y)) 
denote the set of effective variables of /(y). Intuitively, this 
asymmetry allows us more freedom to control the values of 
Pi(x) and 77 (x) such that they can change independently, 
which makes the network behave more like a wireless channel. 
The formal description of the method is presented below: 

Corollary VI.l. Assume all mij{n)'s (i,j = 1, 2, 3) are non- 
zeros, and rjijK.) is not constant. Each unicast session can 
asymptotically achieve one half rate through PBNA if for 
1^1,2, 3, p,(x) ^ 1 and V(?7(x)) ^ V(p,(x)). 

Proof: If the above conditions are satisfied, we must have 
p,(x) ^ e S'. Thus, the theorem holds. 



In Fig. [4a] and 4b we present two examples for which the 
simple method is applicable. As shown in these examples, due 
to edge e, ?7(x) contains effective variables Xcr,je,XeT2^ which 
are absent in the unique form of Pi(x) (i — 1,2,3). Thus, 
by Corollary VI.l each unicast session can asymptotically 



achieve one half rate through PBNA. However, Corollary VI. 1 



doesn't subsume all possible networks for which PBNA is 
feasible. For instance, in Fig. |4c] we show a counter example, 
where V(77(x)) 
applicable. 



V(pi(x)), and thus Corollai-y VI.l 
Nevertheless, it is 



is not 

easy to verify the network 
satisfies the conditions of the Main Theorem, and thus PBNA 
is still feasible. 




Fig. 4. Applying the simple method of Corollary |VI.1| to networks: In the 
top two figures, due to edge e, ??(x) contains coding variables Xa^e, XeT2^ 
which are absent in the unique forms ofpi(x),p2(x) and (x) . Thus, the 
simple method applies to both networks. In contrast, the simple method fails 
in the third network due to V(r;(x)) = V(pi(x)). However, PBNA is still 
feasible for this network according to the Main Theorem. 



2) Reducing S' to S'^: As shown in Section III-D not all 
coupling relations Pi{x) = e S' are realizable due to 

the special properties of transfer functions. Indeed, using the 
two graph-related properties we introduce in Section [V] and 
a recent result of ^1T\ , we can reduce S' to the minimal set 
S', which represents the set of all coupling relations Pi(x) = 



/('?(x)) 



that are realizable. 



9(')(x)) 

The proof consists of three steps. First, we use Linearization 
Property and a simple degree-counting technique to reduce S' 
to the following set 5": 



5r 



ao + ai»y(x) 



e S' : ao,ai,6o,&i e 



(23) 



bo + &i?7(x) 

Note that 5" only includes a finite number of rational func- 
tions. Next, we iterate through all possible configurations of 
ao, ai,bo,bi, and utilize Linearization Property and Square- 
Term Property to further reduce 5" to just four rational 
functions: 

77(x) 



l,r7(x), 1 + ?7(x), 



(24) 



l + 77(x) 

Finally, we use a recent result from [17J to rule out the fourth 
redundant rational function in S2, resulting in the minimal set 
5j' defined in the Main Theorem. The detailed proof is deferred 
to Appendix [B] 

3) Necessity of the Feasibility Conditions: We first show 
that the choices of precoding matrices are not limited to those 
defined in Eq. 7]|9 In fact, given A, B, C, we can always 
construct a precoding matrix Vi such that Eq. (jS]) is satisfied. 
The construction of Vi involves solving a system of linear 
equations defined on ¥g{z): 

r(z)(zC-BA) = (25) 
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i{z)) € F^+^z). Assume 



where r(z) = (ri(z),-- - , 
ro(z) is a non-zero solution to Eq. (25 i. Substitute z with 
7;(x), and we have ?7(x)ro(?7(x))C — ro(r/(x))BA. Finally, 
construct the following precoding matrix 

^ (rj(r;(xi)) r^{vi^')) ■■■ rj(,?(x2"+i))) (26) 

Apparently, Vi satisfies Eq. (j5]l. Hence, each non-zero solution 
to Eq. (25 1 corresponds to a row of Vi satisfying Eq. (j5]l. 
Conversely, it is straightforward to see that each row of Vi 
satisfying Eq. (j5]l corresponds to a solution to Eq. (25 i. 

As an example, consider the case where n — 2 and q — 4. 
Let a be the primitive element of F4 such that a'^ = 1 and 
+ a + 1 = 0. Moreover, let A = I2 and 





/I 






/«2 


a\ 


c = 


a 






1 






w 








a 1 



It's easy to verify that r(z) = {a z + a, z + a, z + az 



satisfies Eq. (25 1. Thus, we substitute z with ri(x^ 
construct — (r^(r/(x^)) r-^(77(x^)) 



and 



r^(ry(x5))). 

Apparently, Eq. (|5]l is satisfied. From this example, we can 
see that given different A,B,C, we can construct different 
precoding matrix Vi, and thus the choices of precoding 
matrices are not limited to those defined in Eq. (|7]l-(|9]). An 
interesting observation is that the above precoding matrix Vi 
is closely related to Eq. (|7]l through a transform equation: 
Vi = V^F, where 



1 


Actually, this observation can be generaUzed to the following 
Lemma. 

Lemma VI.l. Any Vi satisfying Eq. ([5]) is related to 
through the following transform equation 

Vi = GV^F (27) 

where is defined in Eq. F is an (n + 1) x (n + 1) 
matrix, and G is a {2n + 1) x (2n + 1) diagonal matrix, with 
the (i, i) element being fi{ri{x.^)), where fi{z) is an arbitrary 
non-zero rational function in ¥q{z). Moreover, the (n + l)th 
row of FC and the 1st row of FBA are both zero vectors. 

Proof: See Appendix [B] ■ 
Using Lemma |VI. 1| we can prove that when PBNA is 
infeasible via the precoding matrices defined in Eq. (|7])-(j9]) 
due to the presence of a coupling relation, it is infeasible via 
any precoding matrices which satisfy the alignment conditions. 
This implies that the conditions of the Main Theorem are also 
necessary for the feasibility of PBNA. We defer the detailed 
proof to Appendix [B| 



B. rj{-x) Is Constant 

For this case, we use a scheme that is slightly different 



from that of Section III In this scheme, at each sender edge 
ai (i = 1,2,3), we encode source symbol Xi into two encoded 




Fig. 5. An example where r;(x) = 1 and Pi(x) 7^ 1 for i S {1, 2, 3}, and 
thus each unicast session can achieve one half rate in exactly two time slots 
due to Theorem IIV.2I 



symbols, which are then transmitted through the network in 
two time slots. The precoding matrices are defined as follows: 



V2 =Mi3M23l(01 

V3=Mi2M3-i(0i 



(28) 
(29) 
(30) 



where 9i,62 are arbitrary variables. Because T is an identity 
matrix, it can be easily verified that at each receiver, the 
interferences are perfectly aligned, i.e., 

M21V2 = M31V3 
M32V3 = M12V1 
M23V2 = M13V1 

Moreover, the rank conditions can be reformulated to similar 
forms to 



Proof of Theorem IV.2 



First, assume Pi(x) is not 
constant. Using the above construction, the rank condition for 
the ith unicast session is equivalent to 



det 



'^1 P^{^^)el 

92 P»(x2)02 



= ^^102(k(x1)-P,(x2)) ^0 



Thus (i, ^) is feasible through PBNA. Conversely, if Pi(x) 
is constant and thus P; — I2, the rank condition for the ith 
unicast session is violated, and thus ( ^ , ^ , ^ ) is not feasible 
through PBNA. ■ 
In Fig. |5] we show an example of this case. One interesting 
observation is that in this example, each unicast session can 
achieve unit rate by using pure routing. Hence, PBNA doesn't 
provide any advantage over routing in terms of transmission 
rate. 



C. Some mij{:ic) = (i ^ j) 

In this case, since the number of interfering signals is 
reduced, at least one alignment conditions can be removed, 
and thus the restriction on Vi imposed by Eq. (|5]l vanishes. 
Therefore, we can choose Vi freely, and the feasibility 
conditions of PBNA can be greatly simplified. For exam- 
ple, assume 77121 (x) = and all other transfer functions 
are non-zeros. Hence, the alignment condition for the first 
unicast session vanishes. Using a scheme similar to above, 
we set Vi = (6*1 62)^, Y2 = M13M23H6I1 6*2)^ and 
V3 — Mi2M32^(0i 6*2)"^, and thus the interferences at T2 
and T3 are all perfectly aligned. It is easy to see that (i, i, i) 
is feasible through PBNA if and only if pi (x) is not constant 
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(a) ?7(x) = 1 (b) p2(x) = ??(x) 

Fig. 6. Additional examples of coupling relations 



for every i — 1,2,3. Using similar arguments, we can discuss 
other cases. 



VII. Interpretation and Checking of the 
Feasibility Conditions of PBNA 

In this section, we first interpret the coupling conditions as 
defined in the Main Theorem in terms of the network structure. 
Then, based on these interpretation, we present a polynomial- 
time algorithm to check the feasibility of PBNA. We defer all 
the proofs to Appendix |C] 

In the following discussion, we use the following notations: 

• Cij^ki- The number of edges contained in the minimum 
cut which separates {ai,aj} from {rfe,r/}. 



1) Pi(x) = 1 if and only if Ci2,i3 
and only if C 13^12 — 1. 

2) = 1 if and only if Ci2,23 
and only if 023,12 = 1- 

3) P3(x) = 1 if and only if C23,i3 
and only if 6*13^23 = 1- 



= 1; P2(x) = ry(x) if 

= 1; P3(x) = ?7(x) if 



As shown previously in Fig. 3a due to the single bottleneck 
e, the minimum cut separating all senders and receivers equals 
one, and thus Pi(x) — P2(x) = P3(x) = 77(x) = 1. In 
contrast, in Fig. |6bj there is no single bottleneck, but because 
the minimum cut separating 0-2 , from ri , T2 contains only 
one edge e, the coupling relation P2(x) = 77(x) still holds. 

Now we consider the remaining coupling relations, i.e., 
^'i W = 1+!^' ^'2(x) = 1 + 77(x), and p3(x) = 1 + 77(x). It 
turns out that these coupling relations are more complicated 
than the previous ones, and cannot be easily interpreted in 
terms of cut conditions. In fact, each of them is characterized 



by two edges, as specified in Theorem VII. 3 For instance, 
the coupling relation pi (x) = il^^^y^-^ is characterized by two 
edges 0312 and Ck2i3, which form a cut separating ai from ri. 

Theorem VII.3. Assume ri(x) is not constant. The following 
statements hold: 



Cej^ea" The sct of bottlenecks between two edges ei and 

1) Pi(x) 

62. Here, we define an edge as a bottleneck between ei 
and 62 if it forms a cut separating ei from 62. Moreover, 
we use Cij to represent CaiT,- 

aijk'. The last edge of the topological ordering of the 
edges in n Cik- Here, given a a set of edges T, the 
topological ordering of T is a rearrangement of the edges 
of T, such that for any e, e' G T and head{e) = tail{e'), 
e must precede e' in the ordering. 

f3ijk: The first edge of the topological ordering of the 
edges in Cjk n Ca^.^^r^. 

ei||e2: There is no directed path between two edges ei 
and 62. 



l+»)(x) 



if and only if a^u G C12, 0:213 G C13, 
0312||q:213, and {a3i2, 0:213} forms a cut which sepa- 
rates <7i from Ti. 

2) p2(x) = 1 + 77(x) if and only if 0123 G C23, 0:321 e 
C21, oi23||a32i, and {0123,0321} forms a cut which 
separates 02 from T2. 

3) P3(x) = 1 + 77(x) if and only if O231 G C31, O132 G 
C32, 023i||oi32, and {0231,0132} forms a cut which 
separates (73 from T3. 

We have seen an example of coupling relation pi(x) = 



3b 



A. Interpretation of Coupling Relations 

Since the feasibility conditions of PBNA depend on whether 
77 (x) is constant, we first provide interpretation of the coupling 
relation r/(x) = 1, as shown below: 

Theorem VII.l. ^^(x) = 1 if and only if 0213 — 0312 and 

P213 = /3312- 



Consider the example shown in Fig. 3a It is easy to see that 
in this example, 0213 = 0312 = /32i3 = /33i2 = e, and thus 
?7(x) = 1. In Fig. 6a we show another example of 7/(x) = 1. 



In this example, the last edges of C21 H C23 and C31 n C32 
are both ei, and thus 0213 = 0312 = ei. Moreover, the first 
edges of C13 n Ce^Ta and C12 n Ce^ra are both 62, and thus 
1^213 = P312 = 62- Therefore, /^(x) = 1 according to Theorem 

mn 

The coupling relations Pi(x) = 1 and Pi{x) — ri{x) can be 
easily interpreted in terms of the minimum cut between two 
senders and two receivers, as shown below. 

Theorem VII.2. The following statements hold: 



in Fig. 3b As shown in this figure, 62 is the last edge 
of &3 inC32, arraei the last edge of C2inC23- Thus 62 = 0312 
and ei — 0213. Moreover, 62 G C12, ei G C13, ei||e2, and 
{ei, 62} forms a cut separating txi from n. 

B. Checking the Feasibility Conditions of PBNA 

We assume G' is stored as an adjacency list, i.e., for each 
node V G V , we associate it with the set of its incoming edges 
and the set of its outgoing edges. Moreover, we assume all the 
edges in G' have been arranged in topological order. 

The checking process consists of the following steps: 1) 
Check if /^(x) = 1; 2) if 77(x) = 1, check the conditions of 
Theorem |IV.2t 3) otherwise, check the conditions of the Main 
Theorem. In the following discussion, we present the building 
blocks involved in these steps. 

1} Calculating Cee'-' We use Algorithm [T| to calculate the 
set of bottlenecks C^e' which separates e from e'. The algo- 
rithm consists of two steps: 1) Lines 1-3 are used to calculate 
the set of edges traversed by the paths in Pee', denoted by 
Eee' ■ Note that in the reverse BFS algorithm, we start from e' 
and move upwards by following the incoming edges associated 
with each node. 2) Lines 4-11 are used to calculate Cee' ■ In this 
step, we iterate through each edge e G Eee' in the topological 
order. In each iteration, we calculate C, which forms a cut 
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Algorithm 1: Calculate Ce 



1 Use BFS algorithm to calculate the set of edges reachable from 
e, denoted by Ei\ 

2 Use reverse BFS algorithm to calculate the set of edges which 
is connected to e', denoted by E2; 

3 See' -Bi n E2; 

4 Cee' {(T,}, C {cTi}; 

5 for each e G iJee' in the topological order do 

6 C^C-{e}; 

7 for each e' such that head{e) = tail{e') do 

8 I if e' G See' then C ^ C U {e'}; 

9 end 

10 if C contains one edge then Cee' ^ Cee' U C; 



11 end 



Algorithm 2: Check if pi (x) 



l+r,(x) 



1 0312 the last edge of C31 n C32; 

2 0213 the last edge of C21 n C23; 

3 if 0312 ^ C12 or 0:213 ^ C13 then return false; 

4 Use BFS algorithm to check whether 0312 is connected with 
Q213 by a directed path; 

5 if 0312 is connected with C1213 then return false; 

6 Let Gi denote the subgraph of G' induced by 

E' — {Q!312, 0213}; 

7 Use BFS algorithm to check whether ri is connected to cri in 
Gi; 

8 if ri is connected to ai in Gi then return false; 

9 else return true; 



separating e from e'. If C contains only one edge, we then 
incorporate C into Cee'- The running time of the algorithm is 
0{h\E'\), where h is the maximum in-degree of nodes in G". 

2) Checking if r]{x) = I: Using algorithm [T| and Theorem 
we can easily check whether this coupling relation 
holds. First, we calculate C31 n C32, C21 n C23, from which 



VII. 1 



we get the two edges 0312 and a2i3- Then, we calculate 
C12 nCa3j2,^2, Ci3 n Cq2J3, ^3, from which we get ;33i2 and 
^213- Finally, we use Theorem VII. 1 to check if 77(x) = 1 by 
checking whether Q!312 = /3312 and 0213 = /32i3- 

3) Checking if Pi(x) — I or Pi(x) — ?/(x).- Due to 



Theorem VII. 2 we use Ford-Fulkerson Algorithm to check 
these coupling relations. For example, in order to check 
whether pi(x) — 1, we add a super sender node s', which 
is connected to s'^ and s'2 via two directed edges of capacity 
one, and a super receiver node d', to which d[ and are 
connected via two directed edges of capacity one. We then 
use Ford-Fulkerson Algorithm to calculate the maximum flow 
from s' to d', which is identical to 6*12,13. Thus, by checking 
whether C12.13 — 1, we can identify whether pi(x) — 1. 
Similarly, we can check other coupling relations. 

4) Checking if pi{x) = or p2{x) , p^ix.) = l + ?7(x).- 

We use Algorithm |2j to check if pi (x) = j^^^^^^^ . The other 
two coupling relations can be checked similarly. Note that Line 
4 consists of two steps: First, we start from 0:312 and use BFS 
to check if 0213 is reachable from 0:312; then we start from 
0:213 and use BFS to check if 0312 is reachable from 0213. 
The running time of the algorithm is 0{h\E'\). 



VIII. Conclusion and Future Work 

In this paper, we consider the problem of network coding 
for the SISO network with three unicast sessions. We introduce 
how to apply precoding-based interference alignment |14| to 
this network setting. We show that due to network topology, 
there might be dependency relations, called coupling relations, 
between different transfer functions, which can potentially 
make PBNA infeasible. Using two graph-related properties 
and a recent result from [Tfl, we identify the minimal set of 
coupling relations that are realizable in networks. Moreover, 
we show that each of these coupling relations has a unique 
interpretation in terms of network topology. Based on these 
interpretations, we present a polynomial-time algorithm to 
check the feasibility of PBNA. 

Clearly, the scenario considered in this paper (SISO network 
with three unicast sessions) is the simplest meaningful scenario 
in this setting but there still many problems that remain to be 
solved regarding applying interference alignment technques to 
the network setting. One important problem is the complexity 
of PBNA, which arises in two aspects, i.e., precoding matrix 
and field size, and is inherent in the framework of PBNA. 
Another important issue is how much benefits interference 
alignment provides compared with routing schemes. One di- 
rection for future work is to apply other alignment techniques 
with lower complexity to the network setting. The extensions 
to other networks scenarios beyond SISO with more than 
three unicast sessions are highly non-trivial. Other extensions 
of the current framework to more practical scenarios include 
networks with delays and errors. 

Appendix A 
Proofs of Graph-Related Properties 

A. Linearization Property and Square-Term Property 

The following lemma plays an important role in the proof 
of Linearization Property and the interpretation of the coupled 
relations, Pi{x) ~ 1 and Pi{x.) — 77 (x). The basic idea of this 
lemma is that we can multicast two symbols from two senders 
to two receivers via network coding if and only if the minimum 
cut separating the senders from the receivers is greater than 
one. 

Lemma A.l. mQb(x)TOj,g(x) 7^ maq{'x.)mpb{'x.) if and only if 
there is disjoint path pair (Pi, P2) £ Vat x 'Ppq or (P3, P4) € 

Proof We add a super sender s and connect it to s'^ and 
s'p via two edges of unit capacity, and a super receiver d, to 
which we connect d'^^ and d'^ via two edges of unit capacity. 
Thus, the transfer matrix at d is 



M 



TOafc(x) niaqix.) 

mpf,(x) mpq{x) 



It is easy to see dct(M) = mab(x)TOpg(x) — mag(x)TOpb(x). 
Hence, we can multicast two symbols from s to d, i.e., 
det(M) 7^ 0, if and only if the minimum cut separating s 
from d is at least two, or equivalently there is a disjoint path 
pair (Pi,P2) eP 
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(a) 0(62) > 0(63) and (b) 0(62) > 0(63) and (c) 0(62) < 0(63) 
o(ei) < 0(64) o(ei) > 0(64) 

Fig. 7. The construction of H (in the proof of the Linearization 
Property) enabled by Lemma [AJ] (Pi is disjoint with P2) 



The proof of Linearization Property involves finding a sub- 
graph H such that some coding variable appears exclusively in 
the denominator or numerator of /i(x^f ), i.e., /i(xh) restricted 
to H. In fact, due to Lemma |A.1| such subgraph H always 
exists, if ft,(x) is not constant. 

Proof of Linearization Property: In this proof, given 
a path P, let P[e : e'] denote the path segment of P 
between two edges e and e', including e, e'. We arrange 
the edges of G' in topological order, and for e e E' , 
let o(e) denote e's position in this ordering. Moreover, de- 
note ft,i(x) = mQb(x)mpg(x), /i2(x) = maq{:x.)mpi,{x) 

and 



_ /ii(x) 
d(x) 



and (i(x) — gcd(ft.i(x), /i2(x)). Let si(x) 

■^slx) = Hence gcd(si(x), S2(x)) = 1. It follows 

it(x) = csi(x),t;(x) — cs2(x), where c is a non-zero 
constant in F2™ . By Lemma |A.1[ there exists disjoint path 



pair (Pi, P2) e 
consider the first case. 

We arbitrarily select another path pair (P3, P4) G Vaq x Vpb- 
Since Pi, P3 both originate at a a, and P2, P3 both terminate at 
Tg, there exist ei € PiflPg and 62 € P2nP3 such that the path 
segment along P3 between ei and 62 is disjoint with Pi U P2. 
Similarly, there exist 63 e P2 H P4 and 64 G Pi n P4 such that 
the path segment between 63 and 64 along P4 is disjoint with 
Pi U P2. Construct the following two paths: P3' = Pi[o'a : 
ei] U P;^[ei : 62] U P2[e2 : t ] and P{' = P2[(Tp : 63] U P{[e3 : 



64] U Pi [64 : Tb] (see Fig. |7]l. Let H denote the subgraph of 
G" induced by Pi U P2 U P^' U Pi'. 

We then prove that the theorem holds for H. If 0(62) > 
0(63) (Fig. 7a and 7b 1, the variables along P2[e3 : 62] are 
absent in /i2(x/j). We then arbitrarily select a variable x^e' 
from P2[e3 : 62], and write /ii(x/f) as f{x'^)xee' + g(x^), 
where x'^l^ includes all the variables in xh other than x^e'- 
Meanwhile, /i2(x^f) can be written as /i2(x^)- Clearly, Xee' 
will not show up in (i(x^f) and thus it can also be written as 
(i(x'jLf). We then find values for x^, denoted by r, such that 
f{r)h2{r)d{r) ^ 0. Finally, denote cq — cg{r)d^^{r), ci = 
cf{r)d^^{r) and 0-2 = ch2{r)d^^{r) and the theorem holds. 
On the other hand, if 0(62) < 0(63) (see Fig. 7c 1, the variables 
along Pi[ei : 64] are absent in /i2(xij). We then select a 
variable Xee' from Pi[ei : 64]. Similar to above, it's easy to 
see that u(x) and v{x) can be transformed into Cix^e' + cq 
and C2 respectively. 

For the case where (P3, P4) £ Vaq x Vpb is a disjoint path 




Fig. 8. Illustration of Square-Term Property. A term with x^^, introduced 
by (Pi,P2) in the numerator of h{x) equals another term introduced by 
(P3, P4) in the denominator of /i(x). 



pair, we can show that u(x) and w(x) can be transformed into 
C2 and ci^ee' + Co respectively. ■ 
The basic idea of Square-Term Property is to construct a 
one-to-one mapping between the square terms in the numerator 
of h{x.) and those in the denominator of /i(x). 

Proof of Square-Term Property: First, we define two 



sets Qi ^ {(Pi,P2) e Vab X Vp 



tPi(x)tp,(x)} 



and Q2 - {iP3,P4) e Vaq X Vpb : xj^, \ tp,{x)tp,{x)}. 
Consider a path pair (Pi,P2) G Qi- Since the degree of 
Xee' in ipi(x) and tp^lx) is at most one, we must have 
Xee' I ipi(x) and Xee' \ ipal^). Thus e, e' G Pi n Pj. Let 
Pi , Pi be the parts of Pi before e and after e' respectively. 
Similarly, define P2 and P|. Then construct two new paths: 
P3 = Pi^ U {e, e'} U P| and P4 = Pj U {e, e'} U Pf 
(see Fig. |8]l. Clearly, tp^{x)tp^{x) = tp-^{x)tp^{x), 
and thus (P3,P4) G Q2- The above method establishes 
a one-to-one mapping (p : Qi ^ Q2, such that for 

0((Pl,P2)) = (P3,P4), tp,(x)tp,(x) - tp,{x)tp,{x). 

Hence, /i(x) = ^ E(Pi,P2)eQi W^Ps W = 
ii;; E(P3,P4)eQ. Wip4(x) = /2(x). ■ 

B. Other Graph-Related Properties 

In this section, we present other graph-related properties, 
which reveal more microscopic structures of transfer functions. 



and are to be used in the proofs of Theorems VII. 1 and 



VII. 3 Before proceeding, we first extend the concept of 
transfer function to any two edges e, e' G E' , i.e., niee' (x) = 
X^PeP / ^p(x), where Vee' is the set of paths from e to e'. 

The following lemma states that any transfer function 
r7iee'(x) is fully determined by the two edges e, e'. 

Lemma A.2. Consider two transfer functions meje2(x) and 



and 62 — 64. 



I (x) if and only if Ci = 63 



Proof: Apparently, the "if" part holds trivially. Now 
assume ei 7^ 63 or 62 7^ 64. Then, there must be some 
edge which appears exclusively in Veie2 or Ve^e^, implying 
r7iej^e2(x) 7^ TOe3e^(x). Thus, the lemma holds. ■ 
The following result was first proved by Han et al. p7) . 
It states that each transfer function r7iee'(^) can be uniquely 
factorized into a product of irreducible polynomials according 
to the bottlenecks between e and e'. 

Lemma A.3. We arrange the bottlenecks in Cee' in topological 
order: 61,62,- : 6fe> such that e = ei, e' = e/j. Then, 
rnee'(x) can be factorized as mee'(x) = IliLi^ '^eiei_n (x), 
where rrie-ei^ii^) is an irreducible polynomial. 
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In addition, as shown below, any transfer function niee' (x) 
can be partitioned into a summation of products of transfer 
functions according to a cut between e and e'. 

Lemma A.4. Assutne U — {ei,e2,-- - ,ek} is a cut which 
separates e from e! . If ei\\ej for e,; ^ Cj G U, we have 
TOee'(x) = ^JL-^ mee; (x)me;e' (x). Otherwise, the above 
equality doesn 't hold. 

Proof: For E U, let Vl^, denote the set 

of paths in Vee' which pass through e^. Because 
CiWcj for Ci 7^ e U, V\^, is disjoint with V^^^,. 
Hence, mee'(x) = I]i=i Z^Pe-p- , ^p(x)- Note that 

TOee,(x)me.e'(x) = E(Pi .Pa)^^^" xP,^ ^Pl W^Ps W ■ 

Moreover, each monomial ip(x) in mee'(x) corresponds 
to a monomial <Pj(x)tp, (x) in megj (x)TOe.e' (x). Hence, 
TOee; (x)me;e' (x) = X]pep> ip(x), and the lemma 
holds. On the other hand, if some is upstream of e^, 

Ke' n Pie- ^ 0' and thus mee'(x) ^^1 EpeP-^, ^p(x), 

indicating that the lemma doesn't hold. ■ 

Appendix B 
Proofs of Feasibility Conditions of PBNA 

A. Reducing S' to S[ 

In order to utilize the degree-counting technique, we use 
the following lemma. Basically, it allows us to reformulate 
each ^l^l^jj e S' to its unique form such that we can 

compare the degrees of a coding variable in a(x) and /3(x) 
with its degrees in the numerator and denominator of Pi(x) 
respectively. 

Lemma B.l. Let ¥ be a field, z is a variable and y = 
(2/1 ; 2/2 J • • • ,?/fe) is a vector of variables. Consider four non- 
zero polynomials f{z),g{z) € F[z] and s(y),t(y) G F[y], 
such that gcd{f{z), g{z)) — 1 and gcd(s(y), t(y)) = 1. 
Denote d = max{d/,c?g}. Define two polynomials in ¥[y]: 
«(y) = fCl^yiy) and /3(y) = gC^yiy). Then 
gcd(a(y),/3(y)) = l. 

Proof: See Appendix [P] ■ 
We use the following three steps to reduce S' to S'^. 
Step 1: S' 5^'. Assume p,(x) = e 5'. We will 

prove that d = max{df,dg} ~ 1. Let Pi(x) — ?7(x) ~ 
1^ denote the unique forms of Pi(x) and r/(x) respectively. 
Without loss of generality, let f{z) — '^'^^nO.jZ^ , g{z) ~ 
^j^gbjZ^ where akbi ^ 0. We first consider the case where 
I <k and thus d = k. Define the following two polynomials: 

a(x) = /(rKx))i'=(x) = ^'^^a,i'=-^-(x)5^(x) 
/3(x) = <?(77(x))t'=(x) = fe,t'^-^(x)s^(x) 

Due to Lemma [B.1[ we have a{x) = cu(x),/?(x) = cv{x), 
where c in a non-zero constant in F^. Moreover, according 
to Linearization Property, we assign values to x other than a 
coding variable x^e' such that u(x) and w(x) are transformed 



into: 

u{Xee') = CiXee' + Cq v{Xee') = C2 
or u{Xee') = C2 v{Xee') = CiXge' + Cq 

where cq, ci, C2 G and C1C2 7^ 0. We only consider the first 
case. The proof for the other case is similar. In this case, a(x) 
and /3(x) are transformed into a{xee') = cciXee' + ccq and 
P{xee') — CC2 respectively. 

By contradiction, assume d> 2. We first consider the case 
where I < k and thus d = fc. In this case, we have 

Q;(a;ee') = '^'^_^ajt''^^ {Xee')s\Xee') = CCiXge' + CCq 
PiXee') = ^^_,_^bjt''^\Xee')s^ {Xee') = CC2 
Assume s{Xee') = J2'j^oSjXl^, and t{Xee') = I]^=O^J^ee" 

where s^^r' 0. Thus maxjr-, r'} > 1. Note that the degree 
of Xee' in t'^~^ {xee')s\xee') IS fcr' + — r'). We consider 
the following two cases: 

Case I: r 7^ r'. If r > r', da — kr > 2, contradicting that 
do, = 1. Now assume r < r' . Let li and /2 be the minimum 
exponents of z in /(z) and respectively. It follows that 
da = kr' - li{r' - r) = 1 and dp = /cr' - /al?-' - r) = 0. 
Clearly, ?2 > due to d^ = 0. If r > 0, kr' - hir' - r) > 
kr' — l2r' > 0, contradicting d/j = 0. Hence, r = 0, and h — k 
due to = 0. Meanwhile, da = {k—li)r' = 1, which implies 
that li — k — 1 and r' — 1. Thus, z'^^^ is a common divisor 
of f{z) and g{z), contradicting gcd{ f (z), g{z)) — 1. 

Case II: r = r'. Since da = 1 and d^ = 0, all the terms 
in a{xee') and (3{xee') containing x^'g, must be cancelled out, 
implying that 

Y: h,i^^si = t'; Y b, ( j)' = f^g = 

Hence z — |^ is a common divisor of f{z) and g{z), contra- 
dicting gcd(/(z),.g(z)) = 1. 

Therefore, we have proved d — 1 when I < k. Using similar 
technique, we can prove that d = 1 when I > k. This implies 
that can only be of the form . Hence, we 

have reduced S' to 5". 

Step 2: S'l => We consider the coupling relation 
^'i(^) = Wi^- coupling relations p2(x) = 
and p3(x) = 

can be dealt with similarly. Define 

Qifx) = = ™ii(^)"'ii3(x) ^ ^gg^JJJg j}jg characteristic of 

^-^^ ' Pl{x) r?xi3{x)m2i(x) 

Fg is p. Given an integer m, let nip denote the remainder of 
m divided by p. Since 5" only consists of a finite number 
of rational functions, we iterate all possible configurations of 
ao, ai,bo,bi as follows: 

Case I: ^ = where aiao&i^o 0, and aobi ^ 

aibo. For this case, we have pi{xee') = bo+biPi^^^r'^^x"') ■ 
It immediately follows 

, ^ _ apcj - boCoC2 - bQCiC2Xee' 

QlyXee') — 7221/0 1, \ il2 

Oiqx^g, + (2p6iCoCi - aiCiC2)Xee' + OiCq - aiCQC2 
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Let ui{xee'),vi{xee') denote the numerator and denominator 
of the above equation respectively. Assume ui{xee') \ vi{xee') 
and thus x^e' — is a solution to vi{xee') = 0. 

However, vi{^^^^^^) = ^(ao&i - ai^o) ^ 0. Hence, 
ui{xee') t vi{xee')- Thus, by the definition of g'i(x) and 
Square-Term Property, x^^, must appear in ui{xee'), which 
contradicts the formulation of ui{xee')- 

Case II: ^ = °°+°^"^ , where aoOibo 7^ 0. Similar to Case 
I, we can derive 



umn of C and D respectively. Hence, Ci , • • • , c„ are lin- 
early independent and so are di, • • • ,d„. Assume there ex- 
ist fi{z),--- Jn{z) G F2™(z) such that Y."=i fi{z)iz(^t ^ 
di) — 0. Without loss of generahty, assume ji{z) ~ 
for i € {1,2,- •• , n}, where gi{z),h{z) £ F2m[z]. Thus, 
J2"=i9i{z)i^Ci - di) = 0. Let k = maxig{i_2,... ,„}{dgj 



and assume gi{z) — X]f=o 



Then, it follows 



1=1 



gi{z){zci - di) 



^^{ai.iz'+^Ci 



a-LiZ^di) 



aocl 



hcl^le' + (2pfelCoCl - aiCiC2)Xee' + 6icg " aiCoC2 



1=1 



fe-i 



1=0 1=1 



i=l 



a;+i,idi) 



which contradicts Square-Term Property. 



Case III: 



bp 1 
0.1 ')(x) 



= ^^^""^ ^ 0- Thus ^ 

L. Since the coefficient of each monomial in the 



denominators and numerators of pi(x) and /^(x) equals one, 
it follows = ^ = 1. This indicates that ni(x) = -^r^- 

where ao6o^i 0- It follows that 

&OC0C2 - boCiC2Xee' 



Case IV: 



9(^) 
qi{Xee') 



aoc\ 



bic^ + 2pbiCoCiXee' + bic{xl 



Similar to Case I, this also contradicts Square-Term Property. 



Case V: 



where oq 7^ 0. Hence, qi{xe 



contradicting Square-Term Property. 



Case VI: = ag + aiz, where agai 7^ 0. Thus, it follows 



ao 



1, 



Pi(x) = flo + ai?7(x). Similar to Case III, ai 
implying that pi(x) = 1 + ?7(x). 

Case VII: = aiz, where ai =^ 0. Similar to Case III 

9(2) ' ' 

ai = 1 and hence Pi(x) = 77(x) 

can only take the 
Thus, we have 
reduced 5f to S!^. 

Step 3: S2 5-. We note that in Proposition 3 of p7) , 
it was proved that pi(x) 7^ 1 + //(x), p2(x) 7^ ™d 
P3(x) 7^ i-flj^L) ■ Combined with the above results, we have 
reduced ^2 to S'^. 



Therefore, we have proved that 
form of the four rational functions in S2 



In summary, according to Theorem III. 1 if the conditions of 
the Main Theorem are satisfied, the three unicast sessions can 
asymptotically achieve the rate tuple (i, i, i) through PBNA. 

B. Necessity of the Feasibility Conditions 

As shown previously, each row of Vi satisfying the align- 
ment conditions corresponds to a non-zero solution to Eq. 
([25]l. The following lemma reveals that any non-zero solution 



to Eq. (25 1 is linearly dependent on the particular vector 
), which forms each row of the precoding 

matrix V^;. 



Lemma B.2. Eq. ( 25 i has a non-zero solution in F2rn [z\ in 
the form of r(z) ~ (1, z, z^, • • • , z")F, where F is an [n + 
1) X (n + 1) matrix in F2m. Moreover, any solution to Eq. {25) 
is linearly dependent on (1, z, • • • , z")F. 

Proof: Denote D = BA. First, we will prove that 
rank(zC — T)) — n. Let and d^ denote the ith col- 



- E "o.«dj = 

i=l 

Therefore, the following equations must hold: 



ak,iCi = 

i=l 



n 

E"o. 

i=l 



d,; = 



ai+i^^d,) =0 e {0, 



,k-l} 



Thus ai i — for any i e {1, • • • G 
implying fi{z) = 0. Hence, rank(zC — D) = 71. 

Then, there must be an n x n invertible submatrix in 
zC — D. Without loss of generality, assume this submatrix 
consists of the top n rows of zC — D and denote this 
submatrix by E„_|_i. Let b denote the (n + l)th row of 
zC — D. In order to get a non-zero solution to equation 
(25 1, we first fix r„_|_i(z) = —1. Therefore, equation (25 1 
is transformed into (ri(z),--- ,r„(z))E„+i = b. Let 
denote the submatrix acquired by replacing the ith row 
of E„+i with b. Hence, we get a non-zero solution to 
r(z) = (-5^, . . . , -j^, -1). Moreover, f (z) = 



det E„_n) is also a solution. Note that 



^dctE„+i ' ' dctE 

(Set El , • • • , det E„ , - det E„+ 
the degree of z in each detE^ is at most n. Thus, f (z) can be 
formulated as (1, z, • • • , z")F, where F is an (n+ 1) x (n+1) 
matrix. Since rank(zC — D) = n, all the solutions to equation 
( [25| ) form a one-dimensional linear space. Thus, all solutions 
must be linearly dependent on r(z). ■ 
Based on Lemma |B.2[ we can easily derive that each 
Vi satisfying Eq. (|5]l is related to VJ through a transform 



equation, as defined in Lemma VI. 1 



Proof of Lemma VI. 1 ■ Let be the ith row of Vi, which 
satisfies Eq. (|5]). According to Lemma B.2 must have the 
form /i(77(x'))(l, 77(x*), • • • , ?7"(x'))F, where fi{z) is a non- 
zero rational function in F2™(z). Hence, Vi can be written 



as GV^F. Moreover, Eq. ( [25] l can be rewritten as follows: 

(z, z2, • • • , z"+i)FC (1, z, • • • , z")FBA 

The right side of the above equation contains no z"+^, and 
thus the (?i + l)th row of FC must be zero. Similarly, there 
is no constant term on the left side of the above equation, 
implying that the 1st row of FBA is zero. 



_ /(r)(x)) 

holds for the network. We will prove that is violated 



Now assume the coupling relation Pi(x) = e S' 



16 



for n > 1, and thus it is impossible for the three unicast 
sessions to asymptotically achieve the rate tuple (i, i, i) via 
PBNA. Without loss of generality, assume f{z) — X]fe=o '^kz'^ 
and g{z) — J^^Zo^kz'^- Apparently, if rank(Vi) < n + 1, 
3§l is violated. Thus, in the rest of this proof, we assume 
rank(Vi) = ? 



1. By Lemma VI. 1 Vi 



GV^F, 

where F is an (n + 1) x (n + 1) invertible matrix. The 
jth row of Vi is rj = fj{r]{x^)){l,r]{x^),--- , rj" {x.^ ))F . 
Since the (ri + l)th row of FC is zero, we have r^C = 
fj{ri{x^)){l,rj{x^),--- , ?7"^^(x-'))H, where H consists of 
the top n rows of FC and rank(H) — n. Let a = 
(ao,ai, • • • ,a„)^ and b = (6o,6i, • • • ,^'«-i)'^- For i = 1,2, 
we define a' = F^^a and b' = H^^b. It follows 

r,a' = /,(r;(x^))(l,,y(x^),... ,^«(x^))Fa' 
= /,(77(x^))(l,r7(x^),-.. ,ry«(x^))a 
= /,(,7(x^))/(,7(x^-)) 

= p,(x^)/,(ry(x^))(l,,7(x^'), ' ' ' ,^rHx^))b 
= p,(x^")/,(7y(x^))(l,7?(x^"), • • • ,,7"-i(x^))Hb' 
= p,(x^)r,Cb' 

Hence, the columns of (Vi PiViC) are linearly dependent, 
violating IM'^^. Similarly, we can prove the case of i = 3. 

Appendix C 
Proofs of Interpretations of Feasibility 
Conditions of PBNA 

A. 77(x) = 1 

First, note that r;(x) can be rewritten as a ratio of two ratio- 
nal functions 7?(x) = where /„-fe(x) ^ W".'^ W . 

' /312(x)' Ji]K\ I niikix) 

Hence, in order to interpret ?7(x) = 1, we first study the 
properties of /ij7c(x). 

The following lemma is to be used to derive the general 
structure of /yfc(x). Basically, it provides an easy method to 
calculate the greatest common divisor of two transfer functions 
with one common staring edge or ending edge. 

Lemma C.l. The following statements hold: 

1) For 61,62,63 £ E' such that 62,63 are both down- 
stream of 61. Let 6 be the last edge of the topological 
ordering of the edges in Ceie2 nCe^eg- Then meje(x) = 

gcd(TOeie2 (x), m^^e^ (x)). 

2) For 61,62,63 e E' such that 61,62 are both upstream 
of 63. Let 6 be the first edge of the topological or- 
dering of the edges in Ce^e^ H Cgjea- Then rrieesix.) = 



gcd(TOeie3(x),l 



,(x)). 



Proof: First, consider the first statement. By Lemma [A3] 
the following equations hold: TOe^e^lx) = ?Tieie(x)mee2 (x) 
and 



gcd(r 



1-6163 



(x) - 



(x)TOee3(x). ThuS 



.(x) 



,(x),meje3(x)). Assume gcd(mee2 (x), TOgeg (x)) ^ 



1. By Lemma [A.3| there exists bottlenecks 64,65 such that 

me^e^ix) I gcd(mee2(x),mee3(x)). Clcai-ly, 65 e Ceie^nCeieg 

and 65 is downstream of e, which contradicts that e is the last 
edge of the topological ordering of Ce^ea nCe^ea- Hence, we 
have proved that gcd(mee2 (x), meea (x)) = 1, which in turn 





(a) ctijk fiijk (b) aijt = (Sijf, 

Fig. 9. The structure of fijki^) can be classified into two types: 1) aijk ^ 
I3ijk such that fijk{x) is a rational function with non-constant denominator; 
2) ctijk = Pijk such that fijk(.x) is a polynomial. 



implies that meje(x) — gcd(mejg2 (x), Tie^e, (x)). Similarly, 
we can prove the other statement. ■ 
Using the above lemma, fijk{x) can be reformulated as a 
fraction of two coprime polynomials, as shown below. 

Corollary C.l. fijk{x.) can be formulated as 



/ufe(x) 



™<T,,ft,fc(x)m„^^,,^^.(x) 



Too 



(31) 



where gcd(m<,^,0,^., (x)m„^^.,,^^. (x), to„^^,,^,^., (x)) 
Proof: fijk{x) can be calculated as 

"'V.,a.,fc(x)TO„^^,,^^.(x)TOjfe(x) 



/-yfe(x) 



TO^,,a,^.,(x)m„,^.,^^Jx) 
"^a.,7.:r,(x)TOjfc(x) 
™a.,fc,r,.(x) 

m-a,,k,r, (x)m^^,ff,^., (x)m^,^.,,^, (x) 

"^a.,fc,ft,fc(x)TO;3,,fc,r,(x) 

"^<T,,ft,Jx)m„.^.,,r,(x) 



, gcd(TOc,.^,,^,(x),TO„.^,,^^(x)) = 1 and 
,fc,ft,;e(x),TO„^^,,^^(x)) = 1. Meanwhile, 



By Lemma 
thus gcd(mQ, 

gcd(maij,,,/3ij.Jx),TOo.^.,^,.Jx)) = 1. Hence, we must have 
gcd(TO„,,ft,Jx)TO„,^,,^^(x),m„,^.,^/3^^Jx)) = 1. ■ 
According to Corollary C.l the structure of /^^(x) must 



fall into one of the two types, as shown in Fig. [9] In Fig. 9 
<Xijk 7^ Pijk and fijk{x) is a rational function, the denomi- 
nator of which is a non-constant polynomial TOq,.^^^.^^ (x). 
On the other hand, when 0;^^^ S Cjk and thus a^fe = 
(3ijk, as shown in Fig. |9b] fijk{x.) becomes a polynomial 

TO^^.,a„,(x)TO„.^,,^^.(x). 

Moreover, using Corollary |C.1[ we can easily check whether 
two /ij7c(x)'s are equivalent, as shown in the next corollary. 



It is easy to see that Theorem VII. 1 is just a special case of 
this corollary. 

Corollary C.l. Assume i, j, k,i' , k' e {1,2,3} such that i 7^ 
j,j 7^ k and i' ^ j,j ^ k' . /yfe(x) = /i'jfe'(x) if and only if 

ctijk = ai'jk' and jiijk = Pi'jk'- 



Proof: By Corollary C.l if a. 



tijk 



a,,jk' and (3^jk 
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Pi'jk', we must have 
/zi/c(x) = Ujk'{x) 
ctijk = oii'jk' and Pijk 



fijki^) = /i'jfc'(x). Conversely, if 
"ia«,fc,ftj*=W = m„^,^^,^.^,(x). Thus 
/^i'jfc' by Lemma A. 2 ■ 



fi. pi(x) = 1 and pi{x.) = ?7(x) 



Using Lemma A.l we can easily prove Theorem VII.2 
shown below. 



as 



Proof of Theorem | V7/. 2\ Apparently, by Lemma A. 1 and 
the definition of pi(x), pi(x) = 1 if and only if the minimum 
cut separating cri,(T2 from ri and T3 is one, i.e., Ci2,i3 = 
1. In order to interpret pi(x) = ?7(x), we consider gi(x) = 

^ = "u2(x)m3i(x) - "^"'^^ Pi^'') = equivalent to 

(71 (x) = 1. Similarly, using Lemma |aT| it is easy to see that 
Pi(x) = 77(x) if and only if the minimum cut separating CTi, 0-3 
from Ti,T2 is one, i.e., C13 12 — 1- ■ 



C. pi(x) = iq^ly^ and P2{x),p3{x) = I + r]{x.) 

Note that the three coupling relations can be respectively 
reformulated in terms of fijk{x) as follows: 

Wll(x) = /312(x) + /213(x) 
W22(x) = /l23(x) + /32l(x) 
W33(x) = /23l(x) + /l32(x) 

Thus, as shown below, the three coupling relations can also 
be interpreted by using the properties of fijh{x). 

We only prove state- 



Proof of Theorem VII. 3 



ment 1). The other statements can be proved similarly. 
First, we prove the "if" part. Due to q;312 G C12 and 

"213 e Ci3, /312(x) = m„^^a3i2i^)'^a3i2,ri{^) and 
/213(x) = m^i,a2i3(x)TOa2i3,Ti(x). HenCe, /312(x) + /213(x) 

— 771^1 , £1312 (x)^Q3^2 :Ti (x) ^~ ^CTi ,Q213 (x)^a213 ,Ti (x). 

the other hand, because Q3i2||a2i3 and {a3i2,Q;2i3} 
forms a cut which separates ai from ti, mii(x) = 
,0312 (x)m„3j2,ri (x) +TO^i ^a2i3 (x)TOa2i3,Ti (x) by Lemma 



A.4 



Therefore, mii(x) = /3i2(x) + /2i3(x). 
Next we prove the "only if" part. Assume mii(x) = 

/312(x) + /213(x). If ^312 ^ C12 but 0213 ^ C13, /312(x) is a 

rational function whose denominator is a non-constant polyno- 
mial, while /2i3(x) is a polynomial. Hence /3i2(x) + /2i3(x) 
must be a rational function with non-constant denominator, 
and thus mii(x) ^ /3i2(x)-h/2i3(x). Similarly, if q;312 G C12 
but q;213 ^ C13, we can also prove that r7iii(x) 7^ /3i2(x) + 

/213(x). 

Now assume 0312 ^ C12 and 0213 ^ C13. It fol- 

^312 (x)"'q3i2 ■ 



lows that /3i2(x) 

'^<^l,g213(^)'"°213,Tl (^) 



213 V= '-13- 

-^^^ and /2i3(x) = 

'"0312. P312 y^' 

m « "ra ■ Because ?7(x) ^ 1, we have /3i2(x) ^ 

"*a213'P213 ^ J 

/2i3(x), which indicates that 0312 7^ 0213 or /3312 7^ /32i3 by 
Corollai-y C.2 and "1-0312, fei2 (x) 7^ "^a2i3,/32i3 (x)- Therefore, 
by Lemma A. 3 one of the following cases must hold: 1) 
There exists an irreducible polynomial mee'(x) such that 
mee'(x) I ma3i2,/33i2(x) but mee'(x) I TO„3^3,^,^3(x); 2) 
there exists an irreducible polynomial mee'(x) such that 

rnee'(x) \ m„3,2,/33i2(x) but TOee'(x) | 771^213 ,&i3 (x) • 



Consider case 1). Define the following polynomials: 

/(x) = lcm(m„3,,,^3,,(x),ma2i3Ai3(x)) n ^nd /i(x) = 
/(x)/TOa3i2,/33i2(x) and /2(x) = /(x)7m„,^3^^,^3 (x). 
Hence, we have TOee'(x) -f /i(x), TOee'(x) | /2(x), and 

/312(x) + /213(x) = [m^i,/J3i2(x)TOa3,,,^i(x)/i(x) + 

m„i,,32i3(x)ma2i3,Ti(x)/2(x)]//(x). Moreover, due to 
gcd(m„3i2,/33i2(x),TO„i,^3i2(x)ma3,,,^i(x)) = 1, it 
follows that mee'(x) \ mo.1,^312 (x)mQ3j2 (x). This 
implies that mee'(x) | TOo-j^/33j2(x)m„3j2^T-j(x)/i(x) + 
™'Ti,fti3(x)mQ2i3,Ti(x)/2(x). However, mee/(x) | /(x). 
This indicates that /3i2(x) + /2i3(x) is a rational function with 
non-constant denominator Thus r7iii(x) 7^ /3i2(x) + /2i3(x). 
Similarly, for case 2), we can also prove that 

™ll(x) ^ /312(x) + /213(x). 

Thus, we have proved that 0312 G C12 and 
•2213 G C13. It immediately follows that mii(x) — 

™.Ti,a3l2(x)Wa3^2,ri(x) + TO^i^a2l3(x;)»7^a2i3,ri(x). HenCe 

each path P in 'Pctiti either pass through q;312 or 0213, 
implying that {0312, a2i3} forms a cut separating cti from 
Ti. Moreover, according to Lemma A.4 q;312| |q!213. ■ 



Appendix D 

Proofs of Lemmas on Multivariate Polynomials 



In this section, we present the proof of Lemma |B.1| We 
first prove that Lemma B.l holds for the case where s(x) 
and i(x) are both univariate polynomials. In order to extend 
this result to multivariate polynomials, we employ a simple 
idea that each multivariate polynomial can be viewed as 
an equivalent univariate polynomial on a field of rational 
functions. Specifically, we prove that the problem of checking 
if two multivariate polynomials are co-prime is equivalent to 
checking if their equivalent univariate polynomials are co- 
prime. Finally, based on this result, we prove that Lemma 
IB. II also holds for the multivariate case. 



_ A. The Univariate Case 



In the following lemma, we show that Lemma B.l holds 
for the univariate case. 

Lemma D.l. Let ¥ be a field, and z, y are two variables. 
Consider four non-zero polynomials f{z),g{z) e F[z] and 
s{y),t{y) G F[j/], such that gcd(/(z), 5(2)) = 1 and 
gcd(s(y), = 1. Denote d = niax{df , dg}. Define two 

polynomials a{y) = /(fg{)t''(y) and = gCi^)t''iy). 
Then gcd(a(y), /3(y)) = L 

Proof: Assume w{x) = gcd{a{x), I3{x)) is non-trivial. 
Thus we can find an extension field F of F such that 
there exists aig G F which satisfies w{xo) — and hence 
a{xa) = f3{xo) = 0. In the rest of this proof, we restrict 
our discussion in F. Note that gcd{ f (z), g{z)) = 1 and 
gcd{s{x),t{x)) = 1 also hold for F. Assume t{xo) = and 
thus X — xq \ t{x). Since gcd(s(a;), t(a;)) = 1, it follows 
that X — Xq \ s{x) and thus s{xo) 7^ 0. Hence, either 
"^(^0) 7^ or I3{xq) 7^ 0, contradicting that a{xo), P{xo) 

^We use lcm(/(x), (;(x)) to denote the least common multiple of two 
polynomials /(x) and g(x). 
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are both zeros. Hence, we have proved that t{xo) ^ 0. Then 



and 5(fg}: 



we have /(fgl) = 

which implies that z — 
and g{z), contradicting gcd{ f (z), g{z)) 
proved that gcd(a(t/), = 1. 



0, 



is a common divisor of f{z) 
1. Thus, we have 



B. Viewing Multivariate as Univariate 



In order to extend Lemma |D.1| to the muhivariate case, 
we first show that each muhivariate polynomial can be 
viewed as an equivalent univariate polynomial on a field 
of rational functions. Let y = {yi,y2,--' ,yk) be a vec- 
tor of variables. For any i G {1,2,- •• , fc}, define yi = 
{yi, - ■■ , yi-i,yi+i, ■■■ , yk), i.e., the vector consisting of all 
variables in y other than y.^. Note that any polynomial /(y) G 

F[y] can be formulated as /(y) = /o(yO + fi{yi)yi H h 

fp{yi)yf, where each fj{yi) is a polynomial in F[yj]. Because 
F[yi] is a subset of F(yi), /(y) can also be viewed as a 
univariate polynomial in ¥{yi)[yi]. We use f{yi) to denote 
/(y)'s equivalent counterpart in F(y,:)[?/i]. To differentiate 
these two concepts, we reserve the notations, such as "|", 
"gcd" and "1cm" for field F, and append "1" as a subscript to 
these notations to suggest they are specific to field F(y,;). For 
example, for /(y),g(y) G F[y] and u{yi),v{y,) € ¥{y,)[y,], 
5(y) I /(y) means that there exists h{y) G F[y] such that 
/(y) = h{y)g{y), and u{yi) |i v{yi) means that there exists 
■w{yi) e F[yj](yi) such that v{yi) = w{yi)u{yi). 

Lemma D.2. Assume g{yi) S F[yi] and /(y) € F[y] is 
of the form /(y) = Y.%o fj{yr)yl, where fj{yi) G F[y,]. 
Then g{yi) \ /(y) if and only if g{yi) \ fjiVi) for each 
je {0,1,--- 

Proof: Apparently, if g{yi) \ fjijj) for any j G 
{0, 1,-- - g{y.j) I /(y). Now assume ^(yi) | /(y). Thus 
there exists h{y) G F[yJ such that /(y) = g{yi)h{y). Let 
^(y) = TIj=ahj{yi)yl. Hence, it follows that /j(yj) 
hjiyi)giyr) and thus g(yi) | /j(yj). ■ 



The following result follows immediately from Lemma D.2 



Corollary D.l. Let g{yi) and /(y) be defined as Lemma \D.2\ 
Then gcd{g{yi)J{y)) = gcd(5(y,), /o(yO, • • • , /p(yj)). 

Proof: Note that any divisor of g{yi) must be a 
polynomial in F[y,]. Let d{yi) = gcd{g{y,), f{y)) and 
c?'(yi) = gcd(g(yO,/o(yi)r-- ,/p(yi))- By Lemma IdT 



c^(yj) I fjiYi) for any j G {0, 1, • • • ,p}, implying that 
d{yi) I d'{yi). On the other hand, d'{yi) \ /(y), and thus 
d'(yO I d(y,). Hence, d{y,) = d'(y,). ■ 

Corollary D.2. For t e {1,2,- ■■ ,s}, let /t(y) G F[y] he 
defined as ft{y) = J2T=o ftj(yi)yi' where ftj{yi) G F[y,]. 
Let g{y^) G F[yj]. It follows 

gcd(.9(yj),/i(y), • • • ,/t(y)) 
=gcd(5(yj),/io(yi)> • • • ,/ipi(yi),--- , 
/so(yi)> • • • JspAVi)) 



Proof: We have the following equations 

gcd(5(y,),/i(y),--- ,/t(y)) 
=gcd(5(y,), /i(y), • • • ,.9(y»), /t(y)) 
-gcd(gcd(5(y,), /i(y)), • • • ,gcd(.g(y,), /s(y))) 

=gcd(5(y.0,/io(yj),--- ,/ipi(y-*)>-- - , 

=gcd(5(yO,/io(yj),--- ,/ipi(yO,--- , 
/so(yj),--- JspAy^)) 



Lemma D.3. For i G {1,2,--- , s}, Zef at(y),fet(y) G F[y] 
such that bt{y) ^ and gcd(at(y), 6t(y)) = 1. For t G 
{1,2,-- - ,s}, let vt{y) = lcm(6i(y), • • • ,fot(y)). Then we 
have 



gcd ai(y) 



Vs{y) , ~,vs{y) I s 
My)'"-'"^-^^W'"^^^^ 



1 



Proof: We use induction on s to prove this lemma. Appar- 
ently, the lemma holds for s = 1 due to gcd(ai(y), 6i(y)) = 
1. Assume it holds for s — L Thus it follows 



gcd|^ai(y) 
=gcd^ai(y) 
=gcd( ai(y) 



= gcd( ai(y) 
= gcd( ai(y) 



vs{y) 

bi{y)' 

'Vsjy) 

bi{y)' 

'"Ay) 

bi{y)'' 

'"siy) 

bi{y)' 

'"siy) 

bi{y)' 



,as{y)T-^,Vs{y) 
,..(y)^,&.(y)^ 



gcd(a,(y),6^(y)) 
a6-i(y) 
a6-i(y) 
gcd(i>s-i(y) 



^^(y) 
^.(y) 
Vsjy) Vs{y) 
hs-i{yy bs{y) 
Vsjy) 
bs-i{yy 
Vsiy) 



gcd ai(y) 



^'^(y) 
bi{y)-' 



bs{y) 

,as-\{y)i — ^,-ys-i(y) 



=gcd 



vs{y) 
■Vs-\{y) 



gcd ai(y) 



Ws-i(y), 



&.-i(y) 
«^-i(y) 
^^i(y) '■ 
^^^(y) 



as-i(y) 



bs{y) 
«^-i(y) 
hs-x{y) 



bs{y) 



gcd 



Vs-iiy)'' 



(d) 



gcd 



bs{y) 



^gcd(u,_i(y),6s(y)) 
=gcd(l,i;,_i(y)) = 1 



bsiy). 

Vs-i{y), 



ws-i(y) 



gcd(us_i(y),6^(y)) 



In the above equations, (a) is due to gcd(as(y), 6s(y)) 



1; (b) follows from the fact that 



thus 



bs{y) 



Vs-iiy) and 



gcd(i;s„i(y), ]^); (c) follows from the 
inductive assumption; (d) is due to the equality: Vs{y) = 
lcm(.._i(y),6.(y))^ ^„4-|ggf(^„ . ■ 

In general, each polynomial h{yi) G Ffyj)[yj] is of the 
form hiyO = + f^y. + • • • + f^yf, where for 

each j G {0,1,--- ,p}, aj{yi),bj{yi) G F[y,], bj{y,) ^ 0, 
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gcd{aj{yi),bj{yi)) = 1, and ap{yi) ^ 0. Note that for 
each yj which is absent in h{yi), we let aj{yi) = 
and bj{yi) — 1. Moreover, define the following polynomial 
f^hiYi) =lcm(6o(yj),fei(yi),--- ,bp{y,)). 

Corollary D.3. For j e {1, 2, • • • , s}, let fj{y,) e ¥{y,)[y,]. 
Define v{y{) = lcm{nf^{y,)_, ■ ■ ■ , Hf^{yi)) and fj{y) ^ 
'v{yi)fj{yi)- Thus gcd(u(yi),/i(y), • • • ,/s(y)) = 1 

Proof: Assume fj{yi) has the following form: 



Qjo(yi) ^ Qji(y») ^ 

bjoiyi) bji{yi)' 



where for any j E {1,2,-- - ,s} and t E {0,1,--- ,Pj}, 



ajt{yi),bjtiyt 



e F[y,], bjtiy,) + and 



gcd(ajt(yj),&jt(yj)) = 1. Apparently, v{yi) is 
the least common multiple of all bjtiyiYs. Define 
UjtiYi) = ^ "^[yi]- Hence, we have 

/j(y) = Yl't=oajt{yi)ujt{yi)v\- Then it follows 

gcd(z;(y,),/i(y)),--- ,/.(y)) 

= gcd(u(y^),aio(yi)uio(yi)'--- , aipi (yO^ipi (yOr • • , 

aso{yi)usf){yi), ■ ■ ■ ,aspAyi)uspAyi)) 



where (a) is due to Corollary D.2 and (b) follows from Lemma 



Generally, the definitions of division in F[y] and F(yi )[?/;] 
are different. However, the following theorem reveals the two 
definitions are closely related. 

Theorem D.l. Consider two polynomials f{y),g{y) E F[y], 
where g{y) ^ 0. Then g{y) \ f(y) if and only if g{yi) |i /(y^ 
for every i E {1, 2, • • • , k}. 

Proof: The division equation between f{yi) and g{yi) is 
as follows 



fiVi) = h,{yi)g{y,) + n{y,) 



(32) 







where hi{yi),r,{yi) E ¥{y^)[yi], and either ri{yi) = or 
dr- < dg. Due to the uniqueness of Equation (32i, /(y) | ^(y) 
immediately implies that for any i e {1, 2, • • • , k}, ri{yi) 
and thus g(y,) |i f{y,). 

Conversely, assume for every i e {1, • • • , fc}, g{yi) |i f{yi) 
and hence ri{yi) = 0. Denote hi{y) = M/i. (y*)^t(yO- Clearly, 
hi{y) E F[y]. Then, the following equation holds 

M/i.(y»)/(y) = h{y)g{y) 



By Corollary D.3 gcd(/i,j^ (y^), ft.i(y)) = LThus, iXh^y,) 



g{y). Define g{y) 



Define u{y) 



g(y) 



By Lemma 



g(y) 



D.2 



.9(y) e F[y]. 



uiy) 



gcd(/(y),3(y)) e F[y]. It follows that 

ff(y) 



gcd(/(y),5(y)) 

Mfc,(y»)g(y) 

gcd(/i,(y)g(y),/ift^(y,)g(y)) 

Mfe._(y^)g(y) 

ff(y)gcd(/i,(y),/ift^(yi)) 



^ M/».(y»)g(y) 

5(y) 
= /^/i.(yi) 

Note that variable yi is absent in u(y). Because yi can be 
any arbitrary variable in y, it immediately follows that all the 
variables in y must be absent in u{y), implying that u{y) is 
a constant in F. Hence g{y) \ /(y). ■ 
Moreover, in the next theorem, we will prove that checking 
if two multivariate polynomials are co-prime is equivalent to 
checking if their equivalent univariate polynomials are co- 
prime. 

Theorem D.2. Let f{y),g{y) be two non-zero polynomi- 
als in F[y]. Tlien gcd(/(y), ^(y)) — I if and only if 
gcdi(/(y,),g(y,)) = 1 for any i e {1,2, • • • ,k}. 

Proof: First, assume for any i E {1,2,- ■■ ,k}, 
gcdi{ f (yi), g{yi)) — 1. We use contradiction to prove that 
gcd(/(y),g(y)) = 1. Assume u{y) gcd(/(y), g(y)) is 
not constant. Let yi be a variable which is present in u{y). 
By Theorem [m| u{yi) |i f{yi) and u{yi) |i g{yi), which 
contradicts that gcd^{f{yi),g{yi)) = 1. 

Then, assume gcd{f{y),g{y)) = 1. We also use 
contradiction to prove that for any i E {1,2,- ■■ ,k}, 
gcd^{ f (yi), g{yi)) = 1. Assume there exists i E {I,-'' i ^1 
such that v{yi) — gcd f (yi), g{yi)) is non-trivial. Define 
w{y) = liv{yi)v{yi) E F[y]. Clearly, |i and 

w{yi) |i Thus, there exists E V{yi)[yi] such 

that 

fi.Vi) = w{y{)p{yi) g{yi) = w{y,)q{yi) 

Let s{yi) ^ \cm{fip{yi),Hq{y^)). Define p{y) = s{yi)p{y,) 
and q{y) = s{y,)q{y,). Apparently, p{y),q{y) E F[y]. It 
follows that 

s(y»)/(y) = w{y)p{y) s{y,)g{y) = w{y)q{y) 

Then the following equation holds 

s(yi)gcd(/(y),g(y)) = u'(y)gcd(p(y), ^(y)) 

Due to Corollary |d3] gcd(s(yi), gcd(p(y), ^(y))) = 
?cd(s(y,),p(y),9(y)) = 1. Hence s{yi) \ w{y). Let w(y) = 



j^^- According to Lemma 



D.2 



w{y) is a non-trivial poly- 
nomial in F[y]. Thus, w{y)'\ gcd(/(y), ^(y)), contradicting 
icd(/(y),g(y)) = l. ■ 



C. The Multivariate Case 



Now, we are in the place of extending Lemma D.l to the 
multivariate case. 

Note that if we substitute F 



Proof of Lemma B.l 



with F(yi) and gcd with gcdj^ in Lemma |D.1[ the lemma 
also holds. Apparently, f(z),g{z) E F(y,j)[z]. We will prove 
that gcd^{ f {z), g{z)) ~ 1. By contradiction, assume r{z) = 
gcdi(/(z),g(z)) E F(yOW is non-trivial. Let /(z) = ^ 

and g{z) = Clearly, f{z) and g{z) are both non-zero 

polynomials in ]F(yi)[2]. Then we can find an assignment to 
yi, denoted by y*, such that the coefficients of the maximum 
powers of z in r{z), f{z) and g{z) are all non-zeros. Let f(z) 
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denote the univariate polynomial acquired by assigning = 
y* to r(z). Clearly, f(z) is a common divisor of /(z) and g{z) 
in F[z], contradicting gcd{f{z), g{z)) = 1. Moreover, due to 
gcd(s(y),t(y)) = 1 and Theorem[D!2] gcd^{s{yi),t{y,)) = 1. 
Thus, by Lemma D.l gcdj^(Q;(yi), /3(yi)) = 1. Since i can be 
any integer in {1, 2, • • • , k}, it follows that gcd(Q!(y), (3{y)) — 
1 by Theorem |D.2| ■ 
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